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ABSTRACT 

We present a stellar population analysis of the absorption line strength maps for 48 early- 
type galaxies from the SAURON sample. Using the line strength index maps of H/3, Fe5015, 
and Mg b, measured in the Lick/IDS system and spatially binned to a constant signal-to- 
noise, together with predictions from up-to-date stellar population models, we estimate the 
simple stellar population-equivalent (SSP-equivalent) age, metallicity and abundance ratio 
[a/Fe] over a two-dimensional field extending up to approximately one effective radius. A 
discussion of calibrations and differences between model predictions is given. Maps of SSP- 
equivalent age, metallicity and abundance ratio [o/Fe] are presented for each galaxy. We find 
a large range of SSP-equivalent ages in our sample, of which ^40 per cent of the galax- 
ies show signs of a contribution from a young stellar population. The most extreme cases 
of post-starburst galaxies, with SSP-equivalent ages of ^ 3Gyr observed over the full field- 
of-view, and sometimes even showing signs of residual star-formation, are restricted to low 
mass systems (<j e ^ 100 km s _1 or ^2 xl0 10 Af Q ). Spatially restricted cases of young stel- 
lar populations in circumnuclear regions can almost exclusively be linked to the presence of 
star-formation in a thin, dusty disk/ring, also seen in the near-UV or mid-IR on top of an older 
underlying stellar population. 

T he flattened com ponents with disk-like kinematics previously identified in all fast ro- 
tators (Kraino vic et al.b are shown to be connected to regions of distinct stellar populations. 
These range from the young, still star-forming circumnuclear disks and rings with increased 
metallicity preferentially found in intermediate-mass fast rotators, to apparently old structures 
with extended disk-like kinematics, which are observed to have an increased metallicity and 
mildly depressed [a/Fe] ratio compared to the main body of the galaxy. The slow rotators, of- 
ten harboring kinematically decoupled components in their central regions, generally show no 
stellar population signatures over and above the well known metallicity gradients in early-type 
galaxies and are largely consistent with old (^lOGyr) stellar populations. 

Using radially averaged stellar population gradients we find in agreement with Spolaor 
et al. a mass-metallicity gradient relation where low mass fast rotators form a sequence of 
increasing metallicity gradient with increasing mass. For more massive systems (above ~3.5 
xl0 10 M Q ) there is an overall downturn such that metallicity gradients become shallower 
with increased scatter at a given mass leading to the most massive systems being slow -rotators 
with relatively shallow metallicity gradients. The observed shallower metallicity gradients and 
increased scatter could be a consequence of the competition between different star-formation 
and assembly scenarios following a general trend of diminishing gas fractions and more equal 
mass mergers with increasing mass, leading to the most massive systems being devoid of 
ordered motion and signs of recent star-formation. 

Key words: galaxies: bulges - galaxies: elliptical and lenticular, cD - galaxies: evolution 
- galaxies: formation - galaxies: kinematics and dynamics - galaxies: structure - galaxies: 
nuclei 



1 INTRODUCTION 

In hierarchical galaxy formation models with a A-dominated cold 
dark matter cosmology (ACDM), early-type galaxies are predicted 
to have formed through th e merging of smaller galaxies and accre- 
tion events over time (e.g..|Somerville"& Primack 1999l lCole et al.l 
l200d : |Pe Lucia et ai1l2006l : iBower et alj|2006h . Large galaxy sur- 
veys in the local universe as well as at medium redshift have 
shown that galaxies follow a strong bimodality in colour mag- 
nitude space, forming a 'red sequence' and a 'blue cloud' with 
an intermediate region sparsel y populated by galaxies sometimes 
dubbed the 'green va lley' (e.g jBaldrv et alj2004l ; lBeil et alj2004 
Willm eTet alj|2006h . The hierarchical galaxy formation scenario 
has received renewed observational support by the recent find- 
ings of the COMBO-17 and DEEP2 surveys that the stellar mass 



* Email: hkuntsch@eso.org 



of galaxies on the red sequence has increased by a facto r of 2-5 
since redshift z ~ 1 teell et alj|2004l ; iFaber et alj|2007l : but see 
March esini et a i] |2009h . This growth of mass can arise from two 
main sources (1) the growth of objects already on the red-sequence 
(in the case of cold gas accretion they may move from the red se- 
quence into the 'green valley' during the period of star-formation) 
and (2) the addition of galaxies from the blue cloud once star- 
formation is shut off either by consump tion or by the removal of 
the cold gas reservoir dFaber et al. 2007). 

There is observational evidence that the detailed growth of 
mass on the red sequence is a function of luminosity (here used as 
a proxy for total galaxy mass) as well as a function of environment. 
While m ost massive galax i es have already joined the red sequenc e 
at z ~ 1 dTreu et al.ll2005r . lBundv etai]|2006l : Icimatti et al.ll200^) . 
the low mass end of the red sequence seems to be subject to the 
cont inued arrival and re-rejuvenation of galaxies up to the present 
day dTreu et al.l2005MSchawinski et al.l2007l:ISansom & Northeast! 
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20081 : ISmith et al.lboogl ; IShapiro et al1l2010l. but see iTrager et al.1 



2008). Taken together, a scenario emerges where low mass early- 



type galaxies (on the red sequence), irrespective of their environ- 
ment, are predicted to show younger stellar population ages than 
their more massive counterparts. In addition, galaxies residing in 
high-density environments are older and sho w less spread in ste llar 
age than those in less dense regions (see also lTrager et alj |2008). 

Many authors have investigated the scenarios summarized 
above by means of a stellar population analysis of nearby early-type 
galaxies based on optical line strength indices often measured in the 
Lick/IDS system. In most of th e studies galaxies are treated as spa 



tially unresolved sources (e.g., [Bernardi et al 



2006; ISchawinski et alj l200l ITrager et alj 



120061 ; 



2008; 



Yamada ct al 



Thomas et al 



20101) allowing access to either large samples of gala xies, as for 



exam ple in the case of the Sloan Digital Sky Survey dYork et all 
2000), or the analysis of extremely high signal-to-noise ratio (S/N) 
data. Only a limited number of studies hav e investigated line 
strength gradients obtain ed wit h long slits (e.g.jDavies et alj 19931 ; 
ISanchez-Blazquez et alj|2007l ; ISpolaor et al.|l2009l) or radial ly "av- 
eraging the data from integral-field-units dRawle et al.ll2008r) ; spa- 
tially resolved stellar p opulation maps of early-type galaxies re- 
main the exception (e.g., McDermid et alj20 06; Chilingaria rj2009l ; 



IPracv et al.ll2009l;lBlanc et al.ll2009h 

With the advent of high throughput integral field units on 
medium size telescopes, the study of two-dimensional stellar pop- 
ulations maps has become feasible. We have conducted a survey 
of 72 representative early-type (E, SO, and Sa) galaxies using the 
integral-field spectrograph SA URON mounted on the William Her- 
schel Telescope in La Palma teacon et alj|200l1 ; Ide Zeeuw et al.1 
2002). This allowed us to map the stellar and gas kinematics as 
well as a number of stellar absorption line indices up to about one 
effective radius R c for most of the galaxies in the sample. Ana- 
lyzing the global properties of the velocity fields of the 48 E and 
SO galaxies Emsellem et al.1 d2007l . hereafter Paper IX) character- 
ized the population of early-type galaxy with two classes of slow 
and fast rotators, according to their specific (projec ted) angular mo- 
mentu m measured within one effective radius, A^ jKrainovic et all 
(2008, hereafter Paper XII) further analyzed the velocity and ve- 
locity dispersion maps using kinemetry, a generalization of surface 
photometry to th e higher order momen ts of the line-of-sight veloc- 
ity distribution I Krainovic e t al. 2006) and found that 74 per cent 
of Es and 92 per cent of SOs in the SAURON sample exhibit ro- 
tatio n, which ca n be w ell-described by that of a simple, inclined 
disc. IScott et alj J2009h explore the surprisingly tight correlation 
between the local escape velocity V CS c and stellar population pa- 
rameters. 

In this paper, we present for the 48 E and SO galaxies in the 
SAURON sample, estimates of the mean simple stellar popu lation- 
equivalent (hereafter SSP-equivalent; see ITrager et alj2008l) stellar 
population parameters, age, metallicity and abundance ratio [a/Fe] 
for each bin in our maps. The study is based on the line strength 
measurements presented in Ku ntschner et alj (2006J, hereafter Pa- 
per VI), although we provide a short description of an improved 
data-reduction scheme and the resulting updated line strength mea- 
surements in this paper. Compared to many other line strength 
based stellar population studies, our work rests on a limited num- 
ber (3-4) of line strength indices. However, our study has the ben- 
efit of covering a contiguous region of the galaxies, often reaching 
out to approximately one effective radius R c . Furthermore, it com- 
plements and extends the SAURO N survey studies o f the GALEX 
near-UV (hereafter NUV) imaging jjeong et"al . 2009, hereafter Pa- 



per XIII), Spitzer mid-IR imaging jShapiroetal.ll20ld hereafter 
Paper XV) and ionized gas dSarzi et alj2010l hereafter Paper XVI). 

This paper is organized as follows. In Section [2] we present 
the details of our revised data-reduction leading to improved line 
strength measurements. In Section [3] we describe our method to 
estimate stellar population parameters from the available set of 
Lick/IDS indices along with a discussion on the reliability of the 
results. The two-dimensional stellar population maps are presented 
in Section [4] where we first discuss in detail the findings and also 
analyze overall stellar population gradients and scaling relations of 
the aperture measurements with central velocity dispersion and dy- 
namical mass of the galaxies. The conclusions follow in Section[5] 
We present a revised table of circular aperture line strength mea- 
surements within R c /8 and one R c in the Appendix. 



2 REVISED DATA REDUCTION 

The results presented here are based on the observations and data 
reduction described in lEmsellem etafl fe004 hereafter Paper III), 
Paper VI and lSarzi et al.1 yOOg, hereafter Paper V). However, our 
spectral template library which is used to derive the kine matics 
and the gas emission was changed from the Jones library Jjonesl 
1999) to the stars from the recently published MILES stellar li- 
brary ( Sanchez-Blazquez et al. 2006). This new library provides 
better fit s to our spectra, theref ore minimizing template mismatch 
(see also Cappellari et al. 2007). The main consequence of this are 
small, but significant changes to the kinematics; namely, reduced a 
and hi values where the change is most significant in the central 
regions of our sample galaxies and over a larger area for the most 
massive systems. These changes in the kinematics carry through to 
the fully corrected line strength measurements via the internal ve- 
locity broadening corrections and result in slightly reduced values 
for Mg b, Fe5015 and Fe5270g in the affected regions of the maps 
(see Paper VI). 

In an attempt to further improve the emission-line corrections 
to our absorption line strength measurements, we also re-extracted 
the nebular emission fluxes and kinematics using the MILES stel- 
lar library. In the case of emission-line measurements, however, it 
is crucial to provide an underlying fit to the stellar continuum that 
is as physically motivated as possible in order to avoid spurious re- 
sults. For this reason, rather than using a subset of stellar spectra 
from the MILES library for each separate galaxy, we used for the 
entire sample a single temp late library consisting of 48 simple stel- 
lar population models from lVazdekis et alj (2010, based on MILES 
stellar spectra) to which we added a number of stellar templates 
obtained by matching SAURON spectra devoid of emission. Specifi- 
cally, based on the maps of Paper V, we extracted 50 high-S/N spec- 
tra from circular apertures in regions without significant emission. 
We fitted these with the MILES stars over the entire wavelength 
range (i.e. without excluding regions prone to emission), excluding 
only the stellar templates with exceedingly low (< 0.9 A) values 
of the H/3 absorption-line strength. The optimal combination of the 
stellar spectra that best matched each of our emission-free aperture 
spectra constitute each of the empirical templates in our new tem- 
plate library. Care was taken to extract SAURON spectra covering 
the range of absorption-line strengths observed in Paper VI. In this 
respect, our empirical templates approximate as closely as possi- 
ble the spectra of real early-type galaxies, were they unaffected by 
kinematic broadening. 

Although adopting the new template library does not lead to 
a dramatically better sensitivity in the detection of nebular emis- 
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sion, the use of such a physically-motivated set of templates leads 
to more robust emission-line measurements, and therefore a better 
correction of the absorption-line strength indices affected by emis- 
sion (mostly H/3 and Fe5015 in our survey). For further discussion 
of the new emission-line measurements see Paper XVI. 

In summary, the use of the new spectral template library sig- 
nificantly improved the measurements of the line strength indices. 
Main results presented in earlier papers, however, are not affected. 
Specifically, in Paper VI we presented a comparison between the 
isoindex contours of Mg b and the isophotes of the surface bright- 
ness finding that the Mg b contours appear to be flatter for about 
40 per cent of our galaxies without significant dust features. Us- 
ing the revised index measurements, we find that ellipticity mea- 
surements for the Mg b contours have slightly changed, but agree 
within the one a errors with the previ ous measurements. Furt her- 
more, the stellar M/L ratios derived in C appellari et "all teood) re- 
main the same within the quoted errors. 

In Table I A 1 1 of the Appendix we present the updated line 
strength measurements for circular apertures within R c /8 and one 
R c . These values supersede the ones presented in Paper VI. Aper- 
ture line strength measurements changed by less than 17 per cent 
with robustly measured mean offsets of 0.12, —0.37, —0.02, and 
0.00 A for the H/3, Fe5015, Mgb and Fe5270 s indices over the 
(Rc/8) aperture, respectively. For the one R a aperture, the offsets 
are 0.12, -0.31, -0.01 for the H/3, Fe5015 and Mg b indices, re- 
spectively. Note that we also made use of revised estimates of the 
effective radii R c (see Table I A il l which were derived with a R 1 '" 
growth curve analysis from our wide-field MDM (1.3m) imaging 
survey of the SAURON sample (Falcon-Barroso in preparation). 



3 ESTIMATING STELLAR POPULATION 
PARAMETERS 

Estimates of luminosity-weighted stellar population parameters 
such as age, metallicity and abundance ratios can be inferred from 
a set of line strength indices measured on integrated spectra in 
combination with the predictions from stellar population models. 
In this section we begin by presenting a first order overview of 
the ages and metallicities of our sample galaxies using classical 
age/metallicity diagnostic diagrams and then describe the \ 2 tech- 
nique used for the full analysis of our data including the abundance 
ratio variations. 



3.1 A first look with classical techniques 

The classical problem of the degeneracy between age and metal- 
licity changes in an old (Js2 Gyr) integrated stellar population can 
be partially overcome by plotting a combination of m etallicity- 
sensit ive indices against an age-sensitive index (e.g., Worfhey 
1 1994) . This method has been used successfully in the past and we 
will present an overview of our data with this traditional method 
before we employ state-of-the-art techniques for the full analysis 
of our line strengths maps (see Section [3~2l . 

Due to our limited wavelength range (4800 - 5380 A) wc 
can only use H/3 as an age-sensitive index. In order to minimize 
the influence of abundance ratio variations on our first-order age 
estimates derived from the classical index vs index diagrams we 
define composite indices. These indices are constructed such that 
a mean metallicity is measured rather than a metallicity biased to 
one species. The H /3 index shows only a mild dependence on abun- 
dance ratios (e.g., iThomas et al. and can be safely used in 



a first order analysis. We emphasize that these composite indices 
are largely used for presentation purposes when we want to show 
our measurements in the observational plane of index vs index dia- 
grams. 

Within the SAURON wavelength range we measure three dif- 
ferent metallicity indicators: Fe5015, Mg b, an Fe5270s (for details 
see Paper VI). These indices are not pure metallicity indicators, 
but are also sensitive to age changes. Furthermore, non-solar abun- 
dance ratios can severely aff ect the interpretation of these indices if 
not ta ken into account (e.g., iKuntschner et alfeOQll ; IThomas et al.l 
2003). In practice it is sometimes convenient to define metallicity 
indicators which are insensitive to ab undance ratio var iations. The 
first indejQ o f this sort was defined bv lGonzalej jl993l) and further 
improved bv lThomas etail fe003lfl However, the SAURON wave- 
length range does not allow us to measure Fe5335. Therefore, we 
define new, abundance-ratio-insensitive inde x combinations for the 
available indices. We follow the approach bv lThomas et al.l J2003l) 
and also use their models to derive the index definitions. The sub- 
set of the model library we use to define the new indices spans ages 
from 1 to 15 Gyr, metallicities from [Z/H] = —1.35 to 0.35 and 
abundance ratios of [a/Fe] — 0.0 to 0.5. This covers most of the 
range of stellar population parameters which we probe in this sur- 
vey. 

The new abundance ratio insensitive indices are defined as: 



[MgFe52]' 
[MgFe50]' 



0.64 x Mgb + Fe5270 



0.69 x Mgb + Fe5015 



(1) 



(2) 



The scaling factor for the Mg b index was optimized so that the 
mean difference between solar and non-solar models is zero for this 
index. FigureQ] demonstrates that, similar to the [MgFe]' index de- 
fined bv lThomas et aljj2003l) , the new composite indices are rather 
insensitive to abundance ratio changes over the range probed by the 
model predictions. 

In Figure [2] we show the relation between [MgFe52]' and 
[MgFe50]' versus [MgFe]' fo r a sample of well c alibrated Fornax 
early-type galaxies taken from lKuntschneil (2000). There is a good 
linear correlation between the new composite indices and [MgFe]' 
as expected if they are good indicators of mean metallicity. 

Throughout this paper we use index measurements calibrated 
to the Lick/IDS system (see Paper VI and Section |2j. There 
are va rious stellar population models avai l able for this system 
(e.g.. IWorthevl ll994l ; IVazdekis et alJ 1 19971 : IThomas et al] 120031: 
iLee & Worthevl 120051 : ISchiavor] [2007) which agree to first order 
in their predictions of index strength as function of age, metallic- 
ity and abundance ratio, but differ in details. We demonstrate this 
in Figure [3] where we show radial gradients averaged along lines 
of constant surface brightness in the H/3 versus [MgFe50]' plane 
for the 48 early-type galaxies in our sample. The center of each 
galaxy is indi cated by a filled circle. Overplotted are model pre- 
dictions from lSchiavor] j2007lfl for [ a/Fe] = 0.0 on the left and 
model predictions by Tho mas et al] d2003l) for [a/Fe] = 0.0 on 
the right. This diagram provides a convenient overview of our line 



Mgfe x (Fc5270+Fc5335) 



1 [MgFe] 

2 [MgFe]' = ^Mgb x (0.72 x Fe5270 + 0.28 x Fe5335) 

3 We use here the "base" models from Schiavon (2007) with solar abun- 
dance ratios, where the abundance ratio bias of the stellar library has been 
removed. The model predictions were transformed to the Lick system by 
using the offsets described in Table 1 of Schiavon 2007. 



The SAURON project - XVII 5 



o.i 



0.0 



-0.1 



1 [MgFe]' 



o.i 



1 [MgFe52]' 



-o.i 



10 1 



0.0 



-0.1 



o o 



'I I I I I I' 



'I I I I I I' 



[MgFe50]' 



t 



-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 
[a/Fe] 



Figure 1. The dependence of [MgFe]' iThomas et all 120031) . [MgFe52]' 
and [MgFe50]' on abundanc e ratio [q/Fe], Sho wn is the difference be- 
tween the model predictions I Thomas et al. 2003) at solar abundance ratio 
and the model predictions for non-solar abundance ratios. The model li- 
brary used here spans ages from 1 to 15Gyr, metallicities from [Z/H] = 
— 1.35 to 0.35 and abundance ratios of [«/Fe] = 0.0 to 0.5. 



strength measurements while at the same time showing graphically 
the relation between the stellar population parameters, luminosity- 
weighted age and metallicity, as predicted by the models. 

It is evident from the plot and consistent for both stellar pop- 
ulation models that the galaxies span a large range in luminosity 
weighted age with the majority of galaxies clustering at old ages 
(~ 12 Gyr). However, taking a ~9 Gyr isochrone as a dividing line, 
about 40 per cent of our sample galaxies show signs for the pres- 
ence of young stellar populations. Excluding the very youngest ob- 
jects, radial line strength gradients translate into a decline in metal- 
licity with increasing radius. 

Examining the diagram in more detail one can see that most 
galaxies that are dominated by old stellar populations (e.g., age > 
9 Gyr) exhibit roughly constant H/3 absorption strengt h with radius 
(see also Paper VI). Taken at face value, the lThomas et al.l d2003h 




Figure 2. [MgFe52]' an d [MgFe50]' versu s [MgFe]' for a sample of For- 
nax galaxies taken from Kuntschner (2000) are shown as o pen diamonds. 
Overplotted are stellar population models by Thomas et al. (2003) at solar 
abundance ratio. Average error bars for the Fornax data are shown in the 
upper left corner. Note, however, that the errors are correlated due to some 
indices being used in both axis. 



models predict for those galaxies an age gradient such that the cen- 
ter is younger than the outskirts. Moreover, a significant fraction 
of galaxies shows H/3 strengths at large radii whi ch are below the 
model predictions for an age of 15 Gyr. Using the lSchiavonl (2007 ) 
models both effects are mitigated in the sense that there is little or 
no evidence for age gradients for old galaxies and that the majority 
of line strength measurements lie within the range of model predic- 
tions. 

The TMB03 and Sch iavon models bot h use the same index 
sensitivity functions from iKorn et al.l d2005h to correct for abun- 
dance ratio patter n s, how ever, they use different stellar libraries; 
namely the ljonesl dl999l) library in th e Schiavon models and the 
Lick/IDS library JWorthev et all 19941) in TMB03. 

Although the lSchiavonl j2007t) models appear to better fit the 
observational data of the SAURON survey, we emphasize that stellar 
population models are not considered perfect by us and calibration 
uncertainties between the models and observed data remain an is- 
sue (see also Graves & Schiavorj|200ljh. For the remainder of this 
paper we will use the models of Schiavonl d2007h for our stellar 
population analysis. However, we will make specific reference to 
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Figure 3. Radial line strength gradients averaged along isophotes for the 48 early-type galaxies in the SAURON sample. The center of each galaxy is indicated 
by a filled circle and different colours are used for each galaxy. For galaxies with a strong contribution from a young stellar population (H/3 > 2.5 A) we 
give the galaxy name next to the central data-point. W e also indicate t he locati on of NGC 4486 whic h exhibits H/3 absorption strength substantially below the 
models. Overplotted are stellar population models by Schiavon (2007, left) and Thomas et al. 12003, right; hereafter TMB03) for solar abundance ratios. Note 
that the TMB03 models extend to [Z/H] = +0.67, whereas Schiavon models reach only [Z/H] = +0.2. 



iThomas et al.l d2003l) at key points to demonstrate the model inde- 
pendence of our results. 

3.2 The x 2 technique 

In the classical index-index diagrams described above one does 
typically not make use of all available index measurements and 
furthermore it requires an iterative scheme if one is interested in 
abundance ratio estimates. A more elegant way of deter mining stel- 
lar po pulation parameters is the x 2 technique (e.g., IProctor et al.l 
2004). Here the best fitting stellar population model to a given set 
of index measurements is found by minimizing the differences be- 
tween the measured index values and the stellar population model 
predicti ons in a y 2 sense. This approach was already used by our 
team in McDer mid et alj d2006l Paper VIII), and will be employed 
for the remainder of this paper to derive the estimates of age, metal- 
licity and abundance ratio for a given set of observed indices. To 
allow for a reasonable accuracy in the estimates we interpolate 
the Schiavon models to a fine grid providing 340681 individual 
models, spanning [Z/H] = —1.19 to 0.57 (in steps of 0.02 dex), 
age = 0.89 to 17.7 Gyr (logarithmically, in steps of 0.02 dex), and 
[a/Fe] = -0.2 to 0.5 (in steps of 0.01 dexfl The above model 
parameter space covers the majority of our measurements for the 



sample of 48 galaxies. For index measurements which lie outside 
of the stellar population model parameter space, we use the best 
fitting values at the boundary of the interpolated models. However, 
observed index values are generally covered by the model parame- 
ter space and only individual data points may cross the boundaries. 
We estimate confidence lev els on the derived parameters using the 
Ax 2 values (for details see lPress et al.lll992h . 

We use the x 2 technique to derive SSP-equivalent stellar pop- 
ulation parameters for each bin in our maps. We also use this tech- 
nique to derive stellar population estimates over circular apertures 
of one effective radius R a and of radius R c /8, where all the spec- 
tra in these apertures have been summed up and line strength values 
re-derived. For 28/48 galaxies the SAURON data cover less then one 
effective radius R c and we applied aperture corrections to the line 
strength using the average corrections derived in Paper VI where 
necessary (see Table 4 therein). Corrections for the H/3 index are 
negligible and are relatively small for the Fe5015 and Mg b indices 
with maximum corrections of 8 per cent. The fully corrected set of 
indices is given in Table lAll of the Appendix. 

Due to our limited wavelength range, only the H/3, Fe5015, 
and Mg b indices are available for the full field of view. The fourth 
index, Fe5270s, typically covers a smaller fraction of the field-of- 
view (see Paper VI). However, for a circular aperture of 7? c /8 we 



4 The Schiavon (2007) models are provided as function of [Fe/H] and were 
converted by us to [Z/H] as used throughout this paper by a transformation 



kindly provided by R. Schiavon. As a result of the transformation the mod- 
els do not cover the full metallicity range given above at each [a/Fe]. 
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can derive the Fe5270s values for 37 out of the 48 sample galaxies 
(see also Table lATt . This allows us to compare the stellar popu- 
lation parameters derived by using the Fe5015 index or the more 
widely utilized Fe5270 inde>[fl The results of the comparison are 
shown in Figure|4] 

Because the H/3 index is the primary age indicator and is used 
in both derivations, the age estimates show excellent agreement. 
Overall we also find good agreement for the metallicity and abun- 
dance ratio estimates with maximum differences of 0.08 dex and 
0.11 dex for [Z/H] and [a/Fe], respectively. Linear fits taking into 
account errors in both variables are close to the one-to-one rela- 
tion (see dashed lines in Figure However, in order to obtain a 
good agreement for the abundance ratio estimates we used a cus- 
tom built version of the Schiavon models (Schiavon, private com- 
munication) where [Ti/Fe] = 0.5 x [Mg/Fe]. More specifically 
abundance ratios of [C/Fe], [N/Fe], and [Cr/Fe] are fixed to solar 
ratios, whereas [O/Fe], [Mg/Fe], [Ca/Fe], [Na/Fe], [Si/Fe] are var- 
ied with [Ti/Fe] = 0.5 x [Mg/Fe]. Using the standard version 
of the Schiavon models with [Ti/Fe] = [Mg/Fe] together with 
the Fe5015 index leads to somewhat underestimated [a/Fe] ratios 
for the most massive galaxies. From our data alone it is not clear 
whether the Fe5015 index shows a slightly incorrect calibration in 
the fitting functions or if indeed the [Ti/Fe] ratio is not tracking the 
[Mg/Fe] ratio. Further calibrations with e.g., globular cluster data 
are needed to settle this particular issue. Throughout this paper we 
employ the [Ti/Fe] adjusted models in our estimates of the stellar 
population parameters. 

In Figure [5] we carry out a comparison of the stellar popu- 
lation parameters derived wit h our standard s et of indices (H/3, 
Fe501 5 and Mg 6) between the lSchiavonl J2007l) and lThomas et al.l 
(2003) models. In order to cover the parameter space probed in 
this paper we use the line strength measurements of the Re/ 8 
apertures as well as the index maps of a galaxy with old stellar 
populations (NGC 4660) and one where young stellar populations 
dominate (NGC 4150). Overall we find reasonable to good agree- 
me nt between the mod els with the known 'saturation' effects of 
the lThomas et alj d2003l) model at old ages (15 Gyr; see also Fig- 
ure [3}. For estimates of the metallicities [Z/H] and abundance ra- 
tios [a/Fe], we find good agreement between the models over a 
large parameter space. One difference to note though is that for 
the most metal rich populations found in the center of early-type 
galaxies, the Thomas et al. models predict higher metallicities and 
younger ages compared to Schiavon. This effect is mainly caused 
by the different slopes the mod els predict for H/3 as func tion of 
metallicity (see also Figure [3}. iGraves & Schiavonl |2008) carry 
out a mo re d etailed compa rison between the predictions of the 
ISchiavon] and Thomas et alj models using the sample of galaxies 
from Thomas et al. I J2005h and a so find overall agreement, albeit 
finding more significant o ffsets for metallicity and age estimates 
(see also lSmith et alj|2009l) . 

In summary, we conclude that our stellar population estimates 
derived from H/3, Fe5015 and Mgb and the Schiavon models, 
which we will use for all maps in this paper, agree well with the 
estimates derived from the more widely used index combination of 
H/3, Mg b and Fe5270 and the pre dictions from the stellar popula- 
tion models o f Thom as et all J2003ri . 

As final demonstration of the x 2 technique, we show in Fig- 
ure^ how the one R c aperture index values of our sample of 48 



5 We convert our Fe5270g measurement s to Fe5270 by using F e5270 
1.26 X Fe5270 s + 0.06 (for details see lKuntschner et alj2003) . 



early-type galaxies project onto the model grids in planes of [a/Fe] 
vs [Z/H], [Z/H] vs Age and Age vs [a/Fe]. Error ellipses are drawn 
for each galaxy and demonstrate graphically th e well known er- 
ror correlation between [Z/H] and Age (e.g.. iTrager et al.l |2000| : 
iKuntschner et al.ll200ll) . In Table Q] we present our estimates of 
SSP-equivalent age, metallicity and [a/Fe] for circular apertures 
within R c /8 and one R c . 



3.3 Interpretation of derived stellar population parameters 

A word of caution is necessary when making use of the stellar 
population models described above to interpret the observed line 
strength indices. Each of the stellar population models for a given 
age, metallicity and abundance ratio is made of a set of stars that 
all have the same age, metallicity and abundance ratio. Strictly 
speaking, these models can then only be applied to objects, or spe- 
cific regions within objects for which this is also true. Since we 
observe the integrated light along the line-of-sight for our target 
galaxies we clearly violate this condition even if only metallicity 
gradients are considered. The problem is even larger for galaxies 
with signs of multiple star-formation periods. In practice, the de- 
rived stellar parameters, from any of the methods d escribed above, 
are interpreted as SSP-equivalent. ISerra & Tragei] d2007l) investi- 
gated the SSP-equivalent age and chemical composition measured 
from Lick/IDS line strength indices of a large set of two-population 
(old plus young), composite model stellar populations. Their con- 
clusions are that SSP-equivalent estimates are biased such that the 
derived SSP age depends primarily on the age of the young popula- 
tion, while the derived SSP chemical composition is weighted more 
towards that of the old population. Furthermore, in this context, the 
age of a young stellar component and the mass fraction between 
old and young stars are degenerate. 

Changes in observed line strength measurements due to inter- 
stellar extinction, most likely to be present for early-type galax- 
ies with signs of recent star-formation, are negligible. Basic sim- 
ulatio ns show that with a standard extinction curve dFitzpatrickl 
1 19991) and E(B-V) = 0.7, the H/3, Fe5015, Mg6, and Fe5270 in- 
dices change by less than 0.01, 0.05, 0.01, and 0.02 A, respectively, 
compared to dust free models. 

Given the biases and degeneracies described above, we find 
that the broad conclusions on the ages and chemical composition 
of our sample as described here are still valid, but individual de- 
terminations for a given spectrum can be biased in the presence 
of composite stellar populations. Specifically, the SSP-equivalent 
age estimate is strongly biased towards the last star-formation event 
rather than tracing the age of the majority of the stars contributing 
to the integrated spectrum. 



4 RESULTS 

Figures [T0aT -l below present maps of SSP-equivalent age, metallic- 
ity [Z/H] and abundance ratio [a/Fe] of the 48 early-type galaxies 
in the SAURON sample, ordered by increasing NGC number. For 
each galaxy, we show the total intensity reconstructed from the full 
wavelength range of the SAURON spectra (see also Paper III), and 
the two-dimensional distributions of age, [Z/H] and [a/Fe] over- 
plotted with isophotes of the reconstructed image spaced by single 
magnitude steps. We also include the stellar velocity maps previ- 
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Table 1. List of SSP-equivalent age, metallicity [Z/H] and [ct/Fe] estimates within a circular aperture of R c /S and R c . 



Name Age (Gyr) Metallicity [ct/Fe] Age (Gyr) Metallicity [ct/Fe] 

Rc/S Rc/S Rc/S R G Rc Rc 

(1) (2) (3) (4) (5) (6) (7) 



NGC474 


9 


q+0.4 
^-1.2 





.05 ± 0.04 


0.16 ±0.05 


_ „ in 7 
' - A -0.6 


-0.09 ±0.04 


0.17 ±0.04 


NGC 524 


15 


,4-1 r i 
"^-1.4 





03 ± 0.04 


0.20 ±0.04 


_ „4- 1 9 

12.2+jj 


0.01 ±0.04 


0.18 ±0.05 


NGC821 


11 


o+ 1 ■ 1 

•^-1.0 





.07 ±0.02 


0.22 ± 0.04 


11 9 + 


-0.13 ±0.02 


0.25 ±0.04 


NGC 1023 


11 


_l_0_ 5 
'-1.0 


0. 


.09 ± 0.02 


0.19 ±0.04 


1 1 7 +i - ' 


-0.03 ±0.04 


0.26 ±0.04 


NGC 2549 


3 


r*+0.4 
'-0.3 


0. 


.27 ±0.02 


0.15 ± 0.03 


c- 1 -|-u.o 
°- 1 -0.2 


0.03 ± 0.02 


0.16 ± 0.04 


NGC 2685 


5 


1 _|_o.5 
1-0.4 


-0. 


.01 ±0.04 


0.19 ±0.04 


fi 4 +i - u 


-0.25 ±0.04 


0.21 ±0.05 


NGC 2695 


16 


l_n o 

y -1.5 


-0 


01 ±0.04 


0.31 ±0.05 


1 7 7+ u - 1 
A '-'-0.8 


-0.19 ±0.02 


0.37 ±0.06 


NGC 2699 


6 


. in ft 
•^-0.6 





.21 ± 0.02 


0.16 ±0.04 


-4-0 A 

8.5+"-^ 


-0.09 ±0.04 


0.19 ±0.05 


NGC 2768 


8 


„_i_n q 
' y -0.8 





.01 ±0.04 


0.26 ±0.04 


1 1 9 +il1 
AJ " Z -1.0 


-0.19 ±0.02 


0.28 ±0.06 


NGC 2974 


8 


l_n q 

• y -0.4 


0. 


.11 ±0.04 


0.29 ±0.04 


_ „4_0 Q 
Q Q+U-» 

y '°-0.8 


0.01 ±0.02 


0.24 ± 0.04 


NGC 3032 





q-i-u.i 
y -0.1 





.03 ± 0.04 


0.07 ±0.02 


- - -1-0 1 

A - A -o.i 


0.13 ±0.06 


0.13 ±0.02 


NGC 3156 





.-.-Ml 1 

y -0.1 





.05 ±0.02 


-0.01 ±0.02 


~ „_i_n 1 
9 fr 
z ' u -0.1 


-0.35 ±0.06 


0.02 ±0.04 


NGC 3377 


6 


-7+0.7 

• ' -0.6 





.01 ±0.04 


0.24 ± 0.04 


_ 4-0 A 


-0.23 ±0.04 


0.20 ±0.05 


NGC 3379 


14 


'7+0.7 
■ 7 -l'.9 





01 ±0.04 


0.28 ±0.05 


- , _4- 1 A 

1 4 


-0.11 ±0.02 


0.29 ±0.05 


NGC 3384 


6 


_i_q 3 
A -0.8 





.19 ±0.02 


0.13 ±0.04 


_ - -1-0 7 

'■ A -0.6 


-0.01 ±0.04 


0.16 ±0.04 


NGC 3414 


1.3 


, _1_ 1 '2 

'^-1.2 





.03 ± 0.04 


0.34 ± 0.05 




-0.21 ±0.04 


0.28 ±0.06 


NGC 3489 


1 


_|_Q 1 

■ S -0'.2 





.15 ±0.02 


0.10 ±0.02 


„ „_|_0 1 

2-3±g;i 


-0.09 ±0.04 


0.12 ±0.03 


NGC 3608 


11 


■ 7 ±io 





.05 ± 0.04 


0.24 ±0.04 


10.7±°;g 


-0.11 ± 0.02 


0.24 ±0.04 


NGC 4150 


1 


n+0. 1 
'2-0.1 





07 ±0.02 


0.08 ±0.03 


9 1 +U.1 

z - A -o.i 


-0.17 ±0.06 


0.11 ±0.04 


NGC 4262 


11 


1+1.4 


-0. 


.01 ±0.04 


0.33 ±0.05 


- , _4- 1 A 

14.711;* 


-0.19 ±0.04 


0.34 ±0.06 


NGC 4270 


8 


1 +0.4 
■ 1 -1.0 


-0. 


.07 ±0.04 


0.08 ± 0.04 


_|_Q g 

8 - 5 -0^7 


-0.19 ± 0.04 


0.15 ± 0.05 


NGC 4278 


17 


w+O.l 

■'-1.6 


0. 


.01 ±0.02 


0.46 ± 0.05 


-,-.4-14 
A ^ A -1.2 


-0.09 ±0.04 


0.39 ±0.06 


NGC 4374 


11 


_i_ i 4 
'-1.3 


0. 


.01 ±0.04 


0.30 ± 0.04 




-0.15 ±0.04 


0.31 ±0.06 


NGC 4382 


3 


'7+0.7 
■'-0.2 





05 ± 0.04 


0.16 ± 0.04 


- .4-0 ^ 

5.4+"-° 

— 0.5 


-0.11 ± 0.04 


0.21 ± 0.04 


NGC 4387 


10 


o _|_0.5 
• 2 -1.3 


-0 


.03 ± 0.04 


0.20 ±0.04 


- - _4- 1 7 
1 1 7 + ' 


-0.17 ±0.04 


0.26 ±0.06 


NGC 4458 


10 


•'-0.9 


-0 


11 ±0.02 


0.28 ±0.05 


_ _ — 1- 1 1 
1 1 7 + " 


-0.31 ±0.04 


0.34 ±0.06 


NGC 4459 


3 


_+0 2 
• 5 -0^3 


0. 


.17 ±0.02 


0.22 ±0.04 


_ _ _|_ A 7 

7 1 + u - ' 

'• A -0.6 


-0.13 ± 0.02 


0.20 ± 0.04 


NGC 4473 


11 


—+0.6 

'-1.0 





.09 ± 0.02 


0.24 ± 0.04 


- „ „4-1 & 

±z ' z -0.6 


-0.09 ±0.04 


0.27 ±0.04 


NGC 4477 


9 


—+0.9 

•'-0.9 





.03 ± 0.04 


0.24 ± 0.04 


1 1 7 + ' 


-0.13 ±0.02 


0.26 ±0.04 


NGC 4486 


17 


__l_n I 

■ 7 -oi 





.13 ± 0.02 


0.41 ± 0.05 




-0.03 ± 0.02 


0.44 ± 0.05 


NGC 4526 


6 


*-0.8 





19 ±0.02 


0.22 ± 0.04 


lu -'_0.9 


-0.01 ±0.02 


0.28 ±0.05 


NGC 4546 


10 


T+ 1-0 

•'-0.9 





.13 ±0.02 


0.27 ±0.04 


. - _4- 1 1 

1 1 7 ' ' 

1J -'-i.o 


-0.13 ±0.02 


0.31 ±0.05 


NGC 4550 


4 


Q + 0.4 

■^-0.4 


-0 


.09 ± 0.02 


0.21 ± 0.04 


4+0. D 

°'*-0.8 


-0.25 ± 0.04 


0.15 ± 0.05 


NGC 4552 


10 


2+1.0 
' —0.9 
7 +0.B 
'-0.9 





.21 ± 0.04 


0.25 ± 0.04 


12.8+I? 
— l.i 


0.03 ± 0.04 


0.26 ± 0.05 


NGC 4564 


10 


0. 


.21 ±0.02 


0.19 ±0.03 


10 7+ 1 - 6 
iu -'-0.5 


-0.05 ±0.04 


0.26 ±0.04 


NGC 4570 


11 


7+1.1 

'-1.0 





17 ±0.02 


0.17 ±0.03 


14 1+ 1 - 4 


-0.13 ±0.04 


0.27 ±0.04 


NGC 4621 


13 


4 +i.3 





.09 ± 0.04 


0.28 ±0.05 




-0.11 ±0.04 


0.32 ±0.06 


NGC 4660 


12 


9+1.2 

• z -0.6 





.15 ±0.02 


0.24 ± 0.04 


13-411:1 


-0.11 ±0.02 


0.29 ±0.05 


NGC 5198 


12 


9+1.2 
Z -l.l 


0. 


.13 ±0.04 


0.26 ±0.04 




-0.09 ±0.02 


0.29 ±0.04 


NGC 5308 


12 


0+O.6 
• 8 -l.l 





.05 ± 0.04 


0.24 ± 0.04 


13.4111 


-0.07 ±0.04 


0.29 ±0.05 


NGC 5813 


11 


7+1.7 

•'-0.5 





.05 ± 0.04 


0.33 ±0.05 


13.4+°- 7 


-0.11 ±0.04 


0.32 ±0.06 


NGC 5831 


8 


1 +0.4 

' -0.7 


0. 


.11 ±0.02 


0.13 ±0.04 


'-0.7 


-0.11 ±0.04 


0.13 ±0.04 


NGC 5838 


9 


O+0.9 

••^-o^ 





19 ±0.02 


0.20 ±0.03 


11 2+ - 5 


-0.01 ±0.04 


0.27 ±0.04 


NGC 5845 


11 


9+0.5 
-1.0 





.15 ±0.04 


0.20 ±0.04 


11 7+O.6 

1J -'-i.o 


0.05 ±0.02 


0.24 ±0.04 


NGC 5846 


15 


4+0-7 

*-1.4 


0. 


.09 ± 0.02 


0.28 ±0.04 


17 7 +o.i 

1 '-'-1.6 


-0.09 ±0.04 


0.31 ±0.05 


NGC 5982 


8 


O+0.9 

• y -0.8 





.17 ±0.02 


0.15 ±0.03 


10 7 +0 - 5 
±u - ' -0.9 


-0.03 ±0.02 


0.18 ±0.04 


NGC 7332 


2 


o+O.l 

• 8 -o.i 


0. 


.29 ±0.04 


0.10 ±0.03 


7+0.5 
°- ' -0.2 


-0.03 ±0.02 


0.15 ±0.04 


NGC 7457 


2 


o+O.l 
S -0.1 


0. 


01 ±0.04 


0.07 ±0.03 


3 4+ - 2 


-0.13 ±0.04 


0.11 ±0.04 



Notes: (1) NGC number. (2) - (4) SSP-equivalent estimates of age, metallicity [Z/H] and [ct/Fe] within a circular aperture of R c /S. (5) - (7) SSP-equivalent 
estimates of age, metallicity [Z/H] and [ct/Fe] within a circular aperture of one R c . For galaxies with less than one R c coverage we applied the aperture 
corrections given in Paper VI to the line strength indices and then derived the stellar population estimates. 
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Figure 4. Comparison of stellar population parameters derived from R c /S circular aperture line strength measurem ents (see also Table IAU by using the 
indices H/3, Fe5015 and Mg6 or the indices H/3, Mgfe and Fe5270 with the \ 2 method and the [Ti/Fe] adjusted Schiavon 12007) models as described in 
Section l3~2l The axis labels indicate which one of the Fe5270 or F5015 indices is used. The one-to-one relation is indicated by a solid line in each panel 
whereas the dashed line shows linear fits to the data taking into account errors in both variables. 




Schiavon Age (Gyr) Schiavon [Z/H] Schiavon [a/Fe] 

Figure 5. Comparison of stellar population parameters derived from the stellar population models of ISchiavonl fe007h and lThomaset all feOOl labeled 
TMB03) using the indices H/3, Fe5015 and Mgfe. The black diamonds with error bars represent the R c /S circular aperture values of our sample of galaxies. 
The results from the index maps of NGC 4660 and NGC 4150 are shown as red and blue points, respectively. Representative error bars are shown in the lower 
right corner of each panel for typical stellar population values for each galaxy. 



ously shown in Paper 110 in the bottom panel for each galaxy, used 
as reference for our discussion. The maps are all plotted with the 
same spatial scale, and oriented with respect to the SAURON field 
for presentation purposes. The relative directions of north and east 
are indicated by the orientation arrow next to the galaxy title. The 
maximum and minimum of the plotting range is given in the tab 
attached to each parameter map, and the colour bar indicates the 
colour table used. In order to allow for an easy comparison be- 
tween galaxies, the plotting range of the age maps is fixed to the 
range 0.4 to 18 Gyr, where the colour bars reflect a linear scaling. 
The colours are adjusted such that blue to dark green shades cor- 
respond to old stellar populations (10 — 18 Gyr), while younger 
stellar populations are represented by red and yellow shades. Sim- 
ilarly, the plotting range of the [a/Fe] maps is fixed to the range 

6 The kinematics were re-derived using the new stellar template library 
as described in Section [2] However, there are only very minor (At) ^ 
10 km s -1 ) changes in the recession velocities thus the visual impression 
of the maps is identical to Paper III. 



—0.2 to +0.5 for all galaxies. For the [Z/H] maps we use an inde- 
pendent plotting range for each galaxy in order to better visualize 
the metallicity distribution across the field-of-view. 

In Paper VI we found non-physical values of the H/3, Fe5015 
and Mg b line strengths in the outer regions of some galaxies where 
many individual lenslets are averaged to achieve the target S/N 
of 60. Additionally, we removed regions where a simple stellar 
population interpretation is not meaningful; for example the cen- 
tral AGN contaminated region and the jet of NGC 4486 or the re- 
gions affected by neighbour galaxies of NGC 5846 and the regions 
affected by foreground stars. These bins are indicated with grey 
colour in the maps presented in Figures fTOaT -1. 

Whenever we make reference to the dynamical mass of a 
galaxy, we have estimated it with the scalar virial relation 

M d yn = 5.0 x Recrt/G, (3) 

where a c is the luminosity-weighted second velocity moment 
within one effective (half light) radius R c and the factor 5.0 was 
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Figure 6. Examples of A\ 2 space for finding the best-fitting stellar population parameters (SSP-equivalent age, metallicity [Z/H] and [a/Fe]; Schiavon 2007 
models) to the line strengths of all 48 early-type galaxies as measured within a circular aperture of one R c . The best-fitting model is found by minimizing 
the difference between the model predictions and the three Lick/IDS indices H/3, Fe5015 and Mgfe. The three panels show different projections of the three- 
dimensional solution space, with one of the parameters fixed at the best fitting value. The contours indicate the 1<t confidence levels based on Ax 2 for two 
degrees of freedom. The small black dots show all grid points where \ 2 is evaluated. 



calibrated using detailed dynamical models in Paper IV. This mass 
represents 

M dyn « 2 x Mi /2) (4) 

where M1/2 is the total dynamical mass within a sphere containing 
half of the galaxy light. It should not be confused with the much 
larger total galaxy mass within the virial radius. 

4.1 Overview of stellar population maps 

The stellar population maps in Figs llOaH show a wealth of struc- 
tures. Some general trends are apparent and we discuss these in the 
following. 

4.1.1 Age maps 

To aid the discussion in this section and to provide a simple 
overview of the different age map categories, we show in Figure [7] 
radial line strength gradients averaged along isophotes with the pre- 
dictions of stellar population models overplotted. This is similar to 
Figure [3] However, here we have separated the galaxies into four 
groups according to their radial age gradients. For each galaxy we 
indicate by colour whether it belongs to the fast rotators (blue) or 
slow rotators (red; see Paper IX). An overview of the various star- 
formation indicators used in the SAURON survey is given in Ta- 
ble^ where galaxies have been grouped in the same manner as in 
Figure[3] 

About 60 per cent of the 48 galaxies in our sample show 
roughly flat age maps of dark-blue to green colour indicating an 
overall old age (^ 9 Gyr) over the observed field-of-view. Due to 
the linear sampling of the colour map in age, the similar colour 
in our maps hides potentially significant age variations in absolute 
terms (9-18 Gyr). This reflects, however, the increased error in our 
age estimates towards older ages, which is dictated by the age sen- 
sitivity of the HyS index, becoming more uncertain for the same 
index measurement error at older ages (see also Figure [3j- 

In marked contrast, individual galaxies (NGC 3032, 



NGC3156, NGC3489, NGC4150, and NGC7457) show age 
maps indicating a spatially extended, recent star-formation 
episode covering large parts of the observed field-of-view with 
SSP-equivalent ages of less than ~3Gyr (see Figure |7Ji). Our age 
maps show that the youngest stellar populations in these galaxies 
are found in the central parts, with typically steep age gradients. 
Furthermore, ionized gas, optical dust features or even regular 
disks are detected in this group of objects (see Papers V & XVI 
and Table |2j. 

All of the above galaxies are also detecte d in a systematic sin- 
gle dish CO survey of the SAURON sample dCombes et~al]|2007t) 
further supporting the scenario of recent star-formation. Recently, 
the SAURON observations have been extended with auxiliary data 
in the NUV with GALEX imaging (Paper XIII) and in the mid- 
IR with Spitzer IRAC imaging and IRS spectroscopy (Paper XV). 
The NUV - V colour may be used to trace young stars (age < 1 
Gyr) although conservative limits have to be chosen to avoid c on- 
fusion with the UV-upturn phenomenon (e.g., lO'Connelll 19991 Pa- 
per XIII). Similarly, the existence of young stars and even ongoing 
star-formation is closely connected to the presence of polycyclic 
aromatic hydrocarbon (hereafter PAH) emission around 8 /im (Pa- 
per XV). Remarkably, but perhaps not surprisingly, the four galax- 
ies with the youngest SSP-equivalent ages in our sample (age < 
2 Gyr within an aperture of Re/8) clearly exhibiting signs of recent 
star-formation in the NUV and in the mid-IR where data is available 
(see Table[2]l. This suggests residual ongoing star-formation in this 
group of galaxies. Evidence for ongoing star-formation is found in 
NGC 3032 even in our own emission-line analysis of the ionized 
gas (i.e. a very low [OlII]/H/3 ratio) as demonstrated in Paper XVI. 

The next youngest galaxy in our sample NGC 7457 with a cen- 
tral SSP-equivalent age of 2.8 Gyr (R c /8) is only marginally de- 
tected in the UV and not detected in CO or by PAH emission in the 
mid-IR. We interpret the slightly older SSP-equivalent age and the 
absence of PAH emission in the mid-IR as evidence for an evolved 
version of the younger galaxies discussed above. It is interesting to 
note that all of the above mentioned galaxies have velocity disper- 
sions of cr < 100km s" 1 (M dyYL < 1.6 x 1O 1O M ) and thus be- 
long to the low mass end of our sample. Paper XV argues that these 
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kinds of galaxies are the product of recent gas-rich (minor) merg- 
ers where the cessation of star-formation occurs in an 'outside-in' 
manner as the molecular gas is heated or consumed. The roughly 
equal distribution of co- and counter-rotating kinematic structures 
in these systems supports the external origin of the gas (Paper XV). 

Another group of galaxies exhibit a more localized central re- 
gion of young stars embedded within an otherwise older stellar sys- 
tem (see maps and Figure |7p; NGC4459, NGC4382, NGC4526, 
NGC4550, NGC4552, NGC 4477, NGC5838, and NGC524). In 
all but two cases (NGC 4382, NGC 4552), the recent star-formation 
event - with sometimes residual ongoing star-formation seen in the 
NUV, the mid-IR and the SAURON gas emission analysis (see Ta- 
ble |2j - is located in a thin, embedded disc clearly visible in the 
emission line maps and the unsharp-masked optical dust images 
of Paper V. These thin stellar discs are typically co-rotating with 
the main body of the galaxy albeit at the distinctly higher rota- 
tional speeds of disks. The metallicity maps (see also Section l4.1.2t 
of these galaxies show a spatially coincident region of enhanced 
metallicity, indicating the signs of ongoing metal enrichment due 
to disc star-formation. NGC 524 is perhaps the least obvious mem- 
ber of this group of galaxies. It shows only very mild signs of 
younger stars in the centefl but it features a prominent face-on op- 
tical dust disk with a clear spiral pattern (see Paper V). If there are 
young stars associated with this thin disc, the near face-on orien- 
tation may dilute their contribution to the integrated line-of-sight. 
CO is detected in four out of six cases, with recent interferomet- 
ric CO observations of NGC 4459, NGC 4526, NGC 4550 reveal- 
ing molecular gas discs c o-rotating with the stellar discs seen i n 
the SAURON observations dYoung et al.l2008l ; ICr"ocker et alj2009l) . 
Also the Spitzer IRAC observations (Paper XV) clearly support the 
scenario of star-formation in well defined thin disks for this group 
of galaxies with detections for the majority of observed galaxies 
(see Tabled) 

The SSP-equivalent ages derived for the regions affected by 
star-forming disks in the galaxies above mentioned range from 
about 3 Gyr to ages barely below the age of the universe. This, 
of course, is a reflection of the luminosity weighting and involved 
mass fractions. Typically these disks represent only a m inor mass 
fraction of the total galaxy mass (e.g. JCrocker et alj2009h . but con- 
tain contributions o f very young (~ 100-300 Myr) stars (see also 
iKavirai et af]|2009h . We note that these thin, embedded disks with 
recent star-formation typically reside in galaxies of intermediate 
mass range with a c — 160 — 240 km s _1 (M dyn = 1 — 25 x 
10 10 M ; see also Paper XV). 

For the galaxies with circumnuclear star-forming disks or 
rings, Paper XV suggests a scenario where they experience a tran- 
sient period of renewed star-formation on top of an underlying 
much older population. The origin of the gas may be internal but 
also very minor mergers are conceivable. This process is a good 
candidate for producing the central, rapidly rotating and metal rich 
stellar disks seen in many of our fast rotating galaxies (see Pa- 
per XII and Section l4~L2t . 

There are two galaxies with a localized region of young stars 
(NGC 4382 and NGC 4552) lacking a detection of young stars in 
the NUV and the mid-IR. However, optical HST photometry and 
ground based imagin g provide eviden ce for a recent merger in the 
case of NGC 4382 dllauer et alJ kOQS). In NGC 4382, the unusual 



The S/N per observed spatial element was particularly low for this galaxy 
and for radii larger than 20" with large bins the line strength indices are of 
modest quality. 



bow-tie like distribution of the young stars also seen in an enhanced 
region of the metallicity map can be matched to the stellar kinemat- 
ics, which show a marked twist in these areas with a strong depres- 
sion in the velocity dispersion map (see Figure flOfl and Paper III), 
together with a central counter-rotating component as seen in the 
higher spatial resolution data of Paper VIII. We interpret this as a 
kinematically distinct stellar component (see also Paper XII) made 
of stars from a recent, but more than ~ 1 Gyr ago star-formation 
episode. The ratio of mass involved in this star-formation event to 
the total mass of the galaxy may be significantly higher compared 
to the mass ratio of the young stars in the embedded, thin disks seen 
in the other galaxies of this group. 

For NGC 4552 the observed, mild increase in H/3 absorption 
strengths and thus inferred decrease in SSP-equivalent age of the 
central region is not supported by any of the other recent star- 
formation indicators probed by our survey (see Tabled- Although 
this galaxy is known to exhibit a significant UV-upturn (Bureau et 
al. in prep.), it is unlikely that the H/3 increase originates from this 
since none of the other slow rotator galaxies that show a UV-upturn 
(NGC 4374, 4486, 5198. 5846, 5982; Bureau et al. in prep.) exhibit 
a clear central increase in H/3 absorption strength (see Paper VI) 
or corresponding decrease in SSP-equivalent age. Due to the sig- 
nificant amount of gas emission present in the center of this galaxy 
(see Papers V & XVI), it is also possible that the H/3 absorption has 
been somewhat over-corrected. 

Typically star-formation and the resulting younger SSP- 
equivalent ages are found towards the center of early-type galax- 
ies and star-formation takes place in a thin disk/ring morphol- 
ogy (see Figure [7^-c and also Paper XV). However, there is 
also evidence for star-formation outside the nuclear regions in 
two galaxies. (1) NGC 474 shows an age map with a relatively 
old central part (~9Gyr) while the SSP-equivalent age becomes 
younger towards larger radii (r ^ 6"). This observation is sup- 
ported by blue NUV-optical colours at the same radii (see Pa- 
per XIII) and is perhaps connecte d to the famous shell structures 
of NGC 474 dTurnbulletal1ll999h . (2) NGC 2974 shows signifi- 
cant star-formation in a ring like structure, both clearly exhibiting 
signs of recent star-formation in the NUV and in the mid-IR, out- 
side the obs erved field-of-vie w (hereafter FoV) of the SAURON ob- 
servations jjeong et al . 2007, Pap er XV). There is also HI detected 
in a ring like configuration I Wei i mans et al.ll2008l) . We detect only 
mild signs of younger stellar populations. 

Besides the galaxies with strong signs of recent star-formation 
or spatially constrained star-formation connected to thin, embed- 
ded disks, there is a sizable group of galaxies (~ 30 per cent of the 
sample) that show intermediate SSP-equivalent ages (3-11 Gyr) 
over large parts of the observed FoV (in increasing order of central 
Re/8 SSP-equivalent age: NGC 7457, NGC 7332, NGC 2549, 
NGC 2685, NGC 3384, NGC 2699, NGC 3377, NGC 4270, 
NGC5831, NGC2768, NGC5982, NGC4387, NGC4564, 
NGC 5845). With only two ex ceptions (the polar ring galaxy 
NGC 2768 dCrocker et all 120081) and NGC 2685) these galaxies 
are not detected in CO, indicating that they have exhausted their 
molecular gas supply and are now evolving passively. Moreover, 
neither the dust maps nor the NUV or mid-IR indicators show evi- 
dence for recent (< 1 Gyr) star-formation. This is consistent with a 
star-formation episode in these galaxies which is long enough ago 
to be devoid of dust and very young stars, but recent and massive 
enough to be traced by the H/3 index and thus our SSP-equivalent 
age maps. For these galaxies (see also Figure|7};) we find typically 
flat age maps or only mildly increasing SSP-equivalent age as 
function of radius. 
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Table 2. Summary of star-formation indicators in the SAURON survey. 



Name Age {Re/8) a c Flatter Mg b CO NUV mid-IR Optical Ionized Rotator 

(Gyr) (km s ) contours PAH dust gas (fast/slow) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
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Notes: This table provides a summary of the star-formation indicators used in the SAURON survey. We group galaxies in four major categories and sort 
them by increasing SSP-equivalent age in each category. Column 1 lists the galaxy name. Column 2 the SSP-equivalent age over a circular aperture with 
radius R c /8 (see also Table[TJ. Column 3 gives the average velocity dispersion over an aperture of R c . Column 4 provides information on the shape of the 
Mg b contours compared to the isophotes (see Paper VI). Column 5 , 6 and 7 indicate detections for molecular gas (CO), NUV star-formation signatures and 
mid-IR star-formation detections from the studies of Combes et al. (2007), Paper XIII and Paper XV, respectively. The presence of optical dust and ionized 
gas emission is indicated in columns 8 and 9 (see Papers V & XVI). Detections given in brackets ['(Y)'] are marginal, while 'Y (sf)' indicates signs of 
ongoing star-formation as detected by the analysis of the gas emission in Paper XVI. Column 10 gives our classification into slow and fast rotators (see Paper 
IX for details). Since NGC 4550 features two counter-rotating disks it shows little overall rotation and hence is classified as slow rotator. 
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Figure 7. Age/metallicity diagnostic diagrams for our sample of the early-type galaxies in the SAURON survey. For each galaxy a radial gradient averaged 
along isophotes is shown. The center of each galaxy is indicated by a filled circle and open circle for cluster and field gal axies, respectivel y. The colour 
indicates whether a galaxy belongs to the fast rotators (blue) or slow rotators (red). Overplotted are stellar population models bv lSchiavorJb007h . (a) Galaxies 
with an extended, recent star-formation episode; (b) galaxies with signs of a central, spatially constrained, region of younger stellar populations; (c) galaxies 
with milder but spatially extended signs of a younger stellar population; and (d) galaxies consistent with overall old stellar populations. 



Last, but not least we have the group of galaxies that show 
line strength indices consistent with overall old stellar populations, 
at least within the sensitivity of our line strength observations and 
stellar population analysis (Figure |7Jl). Eight out of these twenty 
galaxies are classified as slow rotators in Paper IX and it is in- 
deed the group of galaxies with old stellar populations where we 
find most of the slow rotators (see also Section |4~3l l. Two galax- 
ies (NGC2695 for larger radii and NGC4486 for all radii) show 
H/3 line strength which are below the stellar population model pre- 



dictions and therefore lead to non-physically large SSP-equivalent 
ages. A thorough check of the data-reduction procedure revealed 
nothing special for these galaxies and thus we can offer no explana- 
tion to why the H/3 line strength is so weak for these two galaxies. 

We note that, when explored at higher spatial resolution, 
individual galaxies in the 'old' group can show signs of re- 
cent star-formation. For example, HST imaging of NGC4570 
dvan den Bosch & Emselleml 1 19981) reveals evidence for a bar 
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driven ring of young stars (~2Gyr) at a radius of 1'.'7 at the pre- 
sumed location of a bar inner Lindblad resonance (ILR). 

It is interesting to explore the stellar populations of the 
class of early-type galaxies with kinematically decoupled compo- 
nents (KDCs). We define them here as sub-components that ei- 
ther rotate around a different axis to the main galaxy, or that ro- 
tate around the same axis, but with an opposite sense of rota- 
tion and thus are distinct from the young discs discussed above. 
In Papers III, VIII and XII, 13 galaxies of the SAURON survey 
were identified to harbor KDCs. Broadly speaking there appear 
to be two types of KDCs: kiloparsec-scale structures which are 
found almost exclusively in slow rotating galaxies with predomi- 
nately old stellar populations (NGC3414, NGC3608, NGC4458, 
NGC 5198, NGC 5813, NGC 5831, NGC 5982); and more compact 
KDCs showing a range in SSP-equivalent ages with a preference to- 
wards more recent star-formation episodes (NGC 3032, NGC 4150, 
NGC 4382, NGC 4621, NGC 7332, NGC 7457; see Paper VIII). 
While the central, thin discs discussed above show a prominent 
connection between young stars, optical dust, nebular emission, 
NUV and mid-IR star-formation indicators and enhanced metal- 
licity, it is remarkable that most of the large scale KDCs do not 
show an equivalent, significant signature in the stellar population 
maps. NGC 5831 and NGC 5982 do exhibit mild signs of a con- 
tribution from young stars, but the spatial location of the KDCs 
and that of the young stars is not well correlated. Additionally, the 
metallicity maps (see Section l4.1.2l > also do not show a significant 
correlation with the rotation direction of the KDC, but rather fol- 
low a normal, smooth radial behavior. In summary, any assembly 
and star-formation scenario designed to explain the observed large- 
scale KDCs must account for the absence of significant stellar pop- 
ulation parameter signatures (i.e. age, metallicity or abundance ra- 
tio). 

This can be achieved, for exampl e, if the formation process 
was completed a long time ago (e.g., |PaviesetalJl200lh . How- 
ever, other possibilities are being ex plored too. A recent dynami- 
cal study J van den Bosch et al. 2008) applying triaxial orbit -based 
models to the prominent KDC galaxy NGC 4365 (see also [Bender 
ll988l ; lDavies etalJ200lh shows that the observed KDC is not phys- 
ically distinct from the main body, but rather caused by a superpo- 
sition of orbits. A more systematic study of the dynamics of KDC 
galaxies is in preparation (van den Bosch et al.) but it may well be 
that a significant fraction of the observed large scale KDCs are not 
dynamically distinct subcomponents, but rather are the result of a 
projection of orbits from a more extended structure, thus naturally 
explaining the absence of any distinct stellar population signatures. 

4.1.2 Metallicity and abundance ratio maps 

For galaxies with old stellar populations, the abundance ratio maps 
often show to first order no or only mild positive radial gradients 
with values of [a/Fe] ~ 0.2 — 0.4 (red to yellow colour shad- 
ing). The most significant structures, over and above the generally 
rather noisy appearance of the maps, are typically connected to the 
presence of (very) young stellar populations, which lead to lower, 
more solar values of [a/Fe] (see e.g., NGC3156, NGC4150, and 
NGC 7457). The metallicity maps generally show negative gradi- 
ents with increasing radius, often roughly consistent with the mor- 
phology of the light profile s and consistent with many of the ear- 



lier long slit studies (e.g.. iDavies et alj|l993l ; ICarollo et alii 19931 ; 



ISanchez-Blazquez et alll2007h - Two remarkable outliers from this 
trend exist: NGC 3032 and NGC 4150, showing a central peak in 
metallicity with a surrounding ring of depressed metallicity. These 



two galaxies are among the most extreme examples of a recent star- 
burst in our sample (see TableO. Due to the degeneracies between 
star-formation age and young star mass fractions, as well as the bi- 
ases described in Section [331 the stellar population maps of these 
galaxies are difficult to interpret and a simple interpretation may be 
misleading. 

Mean metallicity and [a/Fe] gradients averaged along 
isophotes are described in more detail in Section |4~2l while we con- 
centrate in this Section on the two-dimensional aspects of the maps. 

In Paper VI we found that the Mg b isoindex contours appear 
to be flatter than the isophotes of the surface brightness for about 
40 per cent of our galaxies without significant dust features (also 
see Table [2]>. This flattening is almost always associated with sig- 
nificant rotation in the galaxies (see also Paper XII) and we inferred 
from this that the fast-rotating component features a higher metal- 
licity and/or an increased [a/Fe] ratio compared to the galaxy as a 
whole. With the stellar populations maps presented in this study we 
confirm that the metallicity is enhanced for fast rotating systems 
such that the iso-metallicity contours appear more flattened than 
the isophotes (prominent examples include NGC 4564, NGC 4473, 
NGC 4621 and NGC 4660; see Figure [§]>. However, the region of 
enhanced metallicity is typically connected with a mild decrease in 
our abundance ratio maps (see also Figure [8landll2t. 

For individual galaxies such as NGC 2549 and NGC 3377 the 
region of increased metallicity is connected to the presence of 
younger stellar populations in a rotating structure with disk-like 
kinematics. For other galaxies there is only mild (e.g., NGC 4564) 
or no (e.g., NGC 4660) evidence for younger populations. It is 
tempting to associate these flattened regions of enhanced metal- 
licity and somewhat lower abundance ratios with the more recent 
star-formation seen in thin disks in galaxies such as NGC 4459 and 
NGC 4526 (see Figure [8]l. Taken together, a scenario emerges in 
which the star-formation event that created the fast rotating disc- 
like structure is distinct from the formation of the galaxy main 
body, which took place at a later stage and on somewhat longer 
timescales, leading to an increased metallicity and reduced abun- 
danc e ratio value . 

lNorrisetal] lf2006) demonstrated in their detailed study of the 
lenticular galaxy NGC 31 15, that a simple spheroid + disc forma- 
tion model, where the old spheroid has [a/Fe]= 0.3 and the some- 
what younger disc shows close to solar abundance ratios, works 
well to explain the observed stellar population trends in this galaxy. 

In the SAURON sample, we find good evidence for this sec- 
ondary star-formation in disk-like structures for low to intermedi- 
ate mass (<r c = 100 — 160 km s" 1 ), fast rotator galaxies such as 
NGC 3377 and NGC 4550 where the young stars are still promi- 
nently seen in the SSP-equivalent age maps affecting large regions 
of the galaxy. Can we also find evidence of this type of secondary 
star-formation in older galaxies? 

For more massive galaxies (<r c > 160 km s" 1 ) the disk 
formation and increased metallicity is typically constrained to a 
central location involving only a minor fraction of the total mass 
of the galaxy (e.g., NGC 4459, NGC 4526; Figure Hi. However, 
examining the [a/Fe] maps more carefully, one can see central 
depressions in about 25 per cent of the galaxies (e.g., NGC 821, 
NGC 1023, NGC 2699, NGC 4270, NGC 4387, NGC 4473, 
NGC 4546, NGC 4564, NGC 4570, NGC 4621, NGC 4660, 
NGC 5838). All of these galaxies are classified as fast rotators 
and exhibit overall old (> 9Gyr) stellar populations, while the 
depression in [a/Fe] is highly correlated with an enhancement of 
the metallicity estimates in the same spatial position. A careful 
examination of the age maps often also reveals a weak signature 
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Figure 8. Prominent examples of fast rotators with a distinct, central disk-like component. For each galaxy we show from the top, the stellar velocity field, the 
Mg b line strength map (see Paper VI), the metallicity map, the abundance ratio map [a/Fe] and the SSP-equivalent age map. The disk-like components can be 
clearly seen in the velocity fields and are reflected in enhanced Mg b and [Z/H] regions and mild depressions in the [a/Fe] maps. NGC 4526 and NGC 4459 are 
examples of currently ongoing residual star-formation in a disk. Bins shown in grey colour for NGC 2549 indicate regions of low S/N ratio and thus doubtful 
line strength measurements. Note the varying spatial scales shown for different galaxies; units are given in arcsec. 



of a contribution from a younger stellar population in these 
central regions as compared to the outer regions. However, the 
most important connection to mention is that all of the galaxies 
with central depressions in [a/Fe] show pinched iso-velocity 
contours towards the center indicating the existence of a central 
fast rotating, disk-like component[f| Remarkably, all galaxies of the 
above group that a re included in the analysis of HST imaging by 
lLauer et alj d2005l . 6/12 galaxies) reveal the signature of a stellar 
disk in the central regions. Figures[9landll2lshow radially averaged 
stellar population gradients for selected galaxies where the central 
depressions in [a/Fe] are more apparent. 

The kinemetry analysis of our velocity maps presented in 
Paper XII suggested that all fast rotators contain flattened, fast- 
rotating components contributing at various levels to the total light 
of the galaxy. For example, there are disk dominated galaxies such 
as NGC 3156 and NGC 2685 and more bulge dominated systems 



8 We note that due to inclination effects, edge-on disks in the center are 
more visible than face-on examples. 



such as NGC 821 and NGC4660 (see Paper XII), the latter showing 
largely old stellar populations. Furthermore, 70 per cent of the sam- 
ple shows evidence for distinct, multiple kinematic components. 

A picture emerges where essentially all fast rotator galax- 
ies have seen a period (often multiple periods) of secondary star- 
formation in kinematically disk-like components on top of an 
older, spheroidal component. This picture is perhaps well known 
and a ccepted for classical SOs with sm all bulge to disk ratios 
(e.g jFisher et alJll99rj : lNorris et af]|2006h . However, we find sub- 
components with disk-like kinematics and distinct stellar popu- 
lations in all fast rotators ranging from galaxies classified as Sa 
to (elongated) ellipticals: (a) central very young, still star-forming 
disks in e.g., NGC 4526 (Sa), NGC 4459 (SO; Papers XIII and XV), 
(b) intermediate age systems devoid of any ongoing star-formation 
but clear detections of a younger, metal rich and [a/Fe] depressed 
component with disk-like kinematics as seen in NGC 2699 (E) and 
NGC 3384 (SB) and (c) in the systems with largely old stellar pop- 
ulations where the star-formation event connected to the disk-like 
kinematics is mostly evidenced by a raised metallicity in combina- 
tion with lower [a/Fe] ratios as seen in e.g., NGC 4621 (E5) and 



1 6 Kuntschner et al. 



NGC4660 (E5). Figure ^presents an overview of selected fast ro- 
tators with metal rich, but [a/Fe] depressed disk-like components 
ranging from still star-forming disks as found in NGC 4526 to over- 
all old stellar populations as seen perhaps in NGC 4660. 

The frequent disk-like structures observed at various scales in 
our sample of fast-rotating early-type galaxies suggests dissipative 
processes play a major role in their formation. The very central 
stellar disks seen in HST imaging and the somewhat larger disk- 
like kinematic structures with increased metallicity and depressed 
[a/Fe] seen in our observations of fast rotators match well with 
the recurring but transient circumnuclear star-formation seen in the 
mid-IR Spitzer observations described in Paper XV. The rotational 
support on larger scales for low mass galaxies, which are also of- 
ten showing signs of younger stellar populations is probably con- 
nected to gas-rich (major) mergers adding a wide spread, fast ro- 
tating younger component to the galaxies (see also [Hoffman et all 
2009). It is not obvious that all fast rotating components on larger 
scales in early-type galaxies can be produced in this way. However, 
recent IFU observations for z — 2 ~ 3 galaxies revealed striking 
evidence for a transient phase of clumpy, tu rbulent disk formation 
fed by rapid cold flows jGenzel et al.ll200^) . The old, large-scale 
disk-like structures seen in our sample of galaxies could well be 
the fossil relic of such a cold flow mode of star-formation seen 
at higher redshift. The rotational support and the relatively short 
timescale 1 Gyr) of the secular evolution in these ancient turbu- 
lent gas-rich disks seem s consistent with our observations (see also 
lElmegreen et al]|2008l) . 

The slow rotators in our sample (excluding NGC 4550 with 
its two counter-rotating disks) mostly show flat [a/Fe] maps with 
metallicity gradients consistent with the isophotes. This supports 
their likely different formation path as already argued in Papers IX, 
X and XII. There are however, four slow rotators which show 
abundance ratio maps indicating a somewhat depressed [a/Fe] 
ratio towards the center: NGC 4552, NGC 5198, NGC 5831 and 
NGC 5982. This group of galaxies, which partially also shows 
weak signs of a contribution from a younger stellar population, may 
indicate that the class of slow rotators is perhaps not as homoge- 
neous in their stellar populations as their kinematics may suggest. 

A few words of caution with regard to the results presented 
so far are due. We note that as shown in Paper VI, the Mg b in- 
dex is a very good observational tracer of the increased metallic- 
ity of the disk-like structures. However, in the presence of young 
stellar populations (SSP-equivalent age ^ 3 Gyr) the Mg6 index 
reacts significantly to age and this causes a decrease in absorption 
strength compared to an old population with the same metallicity. 
For example, NGC 4382 (see Paper VI) shows a central Mg b-index 
depression which renders the Mg b index a less useful tool to de- 
tect metallicity enhancements connected to disk-like kinematics in 
these kind of galaxies. 

Using the Thomas et al. models instead of the Schiavon mod- 
els does not alter our main findings and conclusions discussed 
above. In fact, the signs of distinct stellar populations in the regions 
of disk-like kinematics would be more pronounced in metallicity 
for galaxies with old stellar populations. In Figure [9] we show re- 
sults for radially averaged stellar population estimates derived from 
both stellar population models for NGC 462 1 (for details of the pro- 
cedure see Section l4~2] l. 

The main limitation of our abundance ratio estimates is the 
limited number of indices available for the analysis. New obser- 
vations covering large regions of the galaxies as well as a longer 
wavelength range, thus including more indices, are needed to con- 
firm the above findings. 
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Figure 9. Radial gradients of stellar population estimates for the fast rota- 
tor galaxy NGC 4621. The black dots represent individual bins while the 
solid line connects values averaged along isophotes where radial bins are 
indicated by open diamonds. T he stellar populati on estimates are sho wn as 
derived from Schi avonl d2007l left column) and lThomas et alj b003l right 
column) models. The red dots show the results from high spatial resolution 
OASIS observations (see Paper VIII) where the stellar population estimates 
were derived with iThomas et alj g2003f) models, however, using a larger 
set of Lick/IDS indices, namely H/3, Fe5015, Mgb, Fe5270, Fe5335 and 
Fe5406. The vertical dashed lines indicate a radius of 2". 

4.2 Average stellar population gradients 

Average stellar population gradients can be used to study the 
formation history of early-type galaxies since different forma- 
tion models predict di fferent gradien t s. In a nutshell, mono- 
lithic collapse m odels dCarlberd 1 1984 lArimoto & Yoshiil 1 19871: 
IPipino et al.l2008h predict steep metallicity gradients (ranging from 
S[Z/U]/SlogR = A [Z/H] = -0.35 to -1.0) with metal rich 
centers whereas hierarchical models, following a merger tree, pre- 
dict shallower gradients due to the diluti on of any line strength 
gradi ents existing in the pre-merger units dWhite|[r980l : iKob avashi 
2004). However, adding secondary star-formation, introduced by 
gas rich mergers with a range in degree of dissipation, will typi- 
cally lead to enhanced metallicity over the inner regions and thus 
again steeper metallicity gradients. The metallicity gradients are 
predicted to be strongly correlated with the efficiency of the dis- 
sipative pro cess, and to only we akly depend on the remnant mass 
dKoba vashi 2004; Hop kins et al] 2009a. b). 

In this section we derive robust and simple stellar popula- 
tion gradients for SSP-equivalent age, metallicity and abundance 
ratio by averaging the measurements along lines of constant sur- 
face brightness (isophotes) with equal steps in log flux. This ap- 
proach does not take into account the deviations from e.g., metal- 
licity or [a/Fe] contours with respect to isophotes, but does provide 
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Figure 10a. Maps of the derived stellar populations of the 48 early-type galaxies in the SAURON representativ e sample. The SAURON s pectra have been 
spatially binned to a minimum signal-to-noise of 40 by means of the centroidal Voronoi tessellation algorithm of Cappellari & Copin 1 2003). All maps have 
the same spatial scale and the units of x- and y-axis are arcsec. From top to bottom: i) reconstructed total intensity; ii) SSP-equivalent age distribution (the 
age scale is given in units of Gyr and a linear colour scale); iii) metallicity ([Z/H]) distribution and iv) [aFe] distribution. The bottom panel in each column 
reproduces the velocity maps from Paper III. 
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Figure 10b. 
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Figure 10c. 
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Figure lOd. 
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Figure lOe. 
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Figure lOf. 
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Figure lOg. 
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Figure lOh. 
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Figure lOi. 
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Figure lOj. 
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Figure 10k. 
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Figure 101. 
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valuable estimates of mean gradients which can be investigated 
throughout the whole sample and compared to the predictions of 
simulations. Due to the two-dimensional coverage of our data these 
mean radial gradients can be determined with very good precision. 
The average value in each radial bin is derived after applying a 
3cr clipping algorithm. The radius for each bin is calculated as the 
median major axis radius normalized to the effective radius (see 
Table lAll of the Appendix) along the major axis a c = R c /^/l — e, 
where R c is the effective radius and e is the average ellipticity of 
the galaxy within the SAURON field. 

In Fig. QT]we show the radial profiles of SSP-equivalent age, 
metallicity and abundance ratio as function of \og(R/R c ) as de- 
rived from the data. The following conclusions can be inferred from 
the diagram: (1) If the central parts of a galaxy show old (>9 Gyr) 
SSP-equivalent ages, then age gradients remain typically flat. (2) 
In contrast, the presence of young stellar populations in the centre 
typically leads to age gradients such that older SSP-equivalent ages 
are found at larger radii! This observation supports the scenario of 
potentially multiple, centrally located star-formation events which 
add to the older, underlying stars. (3) Excluding galaxies with sig- 
nificant recent episodes of star-formation (central age ^ 3 Gyr), 
metallicity gradients appear remarkably homogeneous with a mean 
value of -0.28 dex per log(R/ R c ) and a rms of 0.12. There is a 
significant range in central metallicity not obviously connected to 
the age of the stellar population. (4) Abundance ratio gradients are 
typically consistent with a flat distribution or are mildly positive 
gradient (median value of 0.06 dex per \og(R/ R c )). Overall, there 
is a trend in the sense that the oldest galaxies have the strongest 
non-solar abundance ratios (see also Figure(6]l. 

For each galaxy we fit an error-weighted straight line to the 
stellar population values at a given radius such that age gradients 
are defined as 

8 log Age 



A log Age = 



5 log(R/Re 



(5) 



where Age is measured in Gyr. A similar definition is used for 
measurements of the gradients in [Z/H] and [a/Fe], denoted as 
A [Z/H] and A [a/Fe], respectively. Examples of the fits are shown 
in Figure [12] for galaxies representative of the group of fast ro- 
tators with young central disks (NGC4150, NGC4382), "old" 
fast rotators with evidence for an embedded disk-like component 
(NGC 4473, NGC 4564) and galaxies from the group of slow rota- 
tors (NGC4374, NGC 3608). 

The error of the measurements in each surface brightness bin 
is taken as the la scatter of the data-points within this bin after 
applying a 3<r clipping algorithm. We further restrict the fitting 
range to all available data with radii between 2" and R c . The inner 
r < 2" regions are excluded from the fit in order to avoid biases 
due to seeing effects. The overall quality of the fits is good and also 
the assumption of linear gradients per A log 7? is reasonable, with 
the exception of some galaxies with contributions from very young 
stars (see e.g., NGC 4150 and NGC 4382 in Figure[l2j. 

In Figures[T3]and[T4]we explore the connection between stel- 
lar population gradients and central stellar population estimates 
within an aperture of Re/ 8, as well as the connection between stel- 
lar population gradients and global kinematic parameters such as 
mean velocity dispersion within one effective radius cr c , the dynam- 
ical mass Mdyn (see Eq.[3]( and the degree of rotational support as 
measured by \r. 

9 A notable exception to this trend is NGC 474 which shows star-formation 
at larger radii (see also Paper XIII and Section l4.1.U . 
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Figure 11. Overview of radial profiles, averaged along isophotes, of SSP- 
equivalent age, metallicity and [a/Fe] for the 48 early-type galaxies in the 
SAURON sample. The colour coding reflects central R c /8 SSP-equivalent 
age with red referring to old stellar populations and blue to young ones (see 
top plot). 



As was already apparent from FigureQT] old (^ 8 Gyr) galax- 
ies generally do not show significant age gradients (mean gradient 
0.02 ± 0.13; the error reflects the one sigma standard deviation of 
the gradients), whereas young, fast rotator galaxies typically show 
positive age gradients (corresponding to a younger center; mean 
gradient 0.28 ± 0.16). Metallicity gradients are, with two excep- 
tions (NGC 3032, NGC 524), all negative with a median value of 
—0.27. The mean metallicity gradient for old (age ^ 8 Gyr) galax- 
ies is —0.25 ±0.11. There is, however, a trend with \r in the sense 
that fast rotators have on average the steeper metallicity gradients 
with a large scatter at high \r. 

The relation between SAURON line strength index gradients 
as well as stellar popula tion gradients with the local escape veloc- 
ity Vesc was explored in lScott et alj j2009h . A tight correlation be- 
tween Vcsc and the line strength indices is found. For Mg b, this cor- 
relation exists not only between different galaxies but remarkably 
also inside individual galaxies - it is both a local and global correla- 
tion. When index measurements are converted to stellar population 
estimates the galaxies are found to be confined to a good approxi- 
mation to a plane in the four-dimensional parameter space of SSP- 
equivalent age, metallicity ([Z/H]), abundance ratio ([a/Fe]) and 
V C sc- The observed tight connection between stellar populations 
and the gravitational potential, both locally and globally, presents a 
strong constraint on galaxy formation scenarios. 
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Figure 12. Examples of radial profiles, averaged along isophotes, of SSP-equivalent age, metallicity and [«/Fe] estimates (open diamonds). We fit an error- 
weighted straight line (orange solid line) to each profile for all data available between 2" (vertical dashed lines) and R c . See text for details. Small black points 
represent individual bins in the stellar population maps. The bottom panels show the kl parameter - a measure of the maximum rotation velocity - from the 
kinemetry analysis of Paper XII. 




[Z/H] (R./8) Age (R e /8) [Gyr] [«/Fe] (R e /8) 

Figure 13. Age gradient, metallicity gradient and [a/Fe] gradient versus metallicity, age and [a/Fe] measured within a circular aperture of R c /S. Fast and 
slow rotators are indicated by blue triangles and red circles, respectively. NGC4550, classified as slow rotator due to two counter-rotating disks, is shown as 
orange filled s ymbol. Cluster and field galaxies are shown as filled and open symbols. In the top middle panel we show in grey the region occupied by the 
simulations of Hopkins et al. (2009a, their Figure 29). 
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Figure 14. Age gradient, metallicity gradient and [a/Fe] gradient versus logcr c , dynamical mass M^ yn an d logA^ for the 48 early-type galaxies in the 
SAURON sample. Fast and slow rotators are indicated by blue triangles and red circles, respectively. Cluster and field galaxies are shown as filled and open 
symbols. NGC4550, classified as slow rotator due to two counter-rotating disks, is shown as orange filled symbol. The one slow rotator at small log <x e is 
NGC 4458. NGC 474, which shows star-formation at larger radii and thus shows a significantly negative age gradient (see also Paper XIII and Section l4. 1 . U is 
labeled in the Age gradient panels. There is a break in the metallicity gradient vs M,j yn relation at a mass of ~ 3.5 X 1O 1O M first detected bv lSpolaor et alj 
120091) and in dicated by the thick grey line. The grey shaded area in the left middle panel shows the approximate region occupied by the simulations of 
iHopkinset all (2009a, their Figure 29). 



ISpolaor et ID (2009) noticed a mass-metallicity gradient 
relation for early-type galaxies where low mass galaxies 
(<C3.5x 1O 1O M0) form a tight relation of increasing metallicity 
gradient with increasing mass. For galaxies above the transition 
mass there is a clear downturn in metallicity gradient with increased 
scatter. Our analysis is consistent with these observa tions (see Fig- 
ure[l4] middle panel, but also see lRawle et alj|201ol) . The fraction 
of slow rotators increases with increasing mass and the most mas- 
sive galaxies are slow rotators (see Paper IX). Furthermore, there is 
a weak trend among the slow rotators in the sense that the galaxies 
with the smallest A_r exhibit the shallowest metallicity gradients. 
This trend clearly needs confirmation with a larger sample of slow 
rotators. 

These observations suggest that the formation and subsequent 
evolution of low-mass early-type galaxies (which are dominated by 
fast rotators) leads to a connection between galaxy mass and metal 
enrichment such that the steepest metallicity gradients are achieved 
for masses of ~3.5xlO lo M0. All these galaxies are rotationally 
supported and feature strong disk-components suggesting gas-rich 
mergers and accretion to play an important role. We emphasize that 
low-mass galaxies with young and old SSP-equivalent ages fol- 
low this relation and thus secondary star-formation due to gas-rich 
(minor) mergers (see Section 14.1.11 Paper XV) does not destroy 
the trend. For more massive galaxies (above ~3.5xlO lo M0) the 
relative influence of the various star-formation and assembly pro- 
cesses which determine the galaxies' properties seem to change. 
We observe a large scatter of metallicity gradients with an overall 



downturn with mass and increased fraction of slow rotators towards 
the most massive systems (> 1O 11 M0). The observed scatter may 
be the result of competing processes such as minor gas-rich merg- 
ers building disk-like structures over extended regions with rather 
shallow metallicity gradients, circumnuclear star-formation in rings 
and disks leading to steeper metallicity gradients, particularly in the 
center, and gas-poor (major) merging removing the signature of or- 
dered motion and also weakening metallicity gradients in the most 
massive remnants. Consistent with this picture, we find the most 
massive systems to be slow rotators that exhibit relatively shallow 
metallicity gradients and are devoid of any signs of recent star- 
formation (see also Paper XV). Brightest cluster galaxies, which 
are u nderrepresented in our sample, may d eviate from the trend 
(e.g.. iBrough et alj|2007l : ISpolaor et alj|2009l) . However, the anal- 
ysis of the projected axial ratio distribution of quiescent galaxies 
in the Sloan Digital Sky Survey ( Ivan der Wei et alj|2009h demon- 
strates that galaxies above a mass of ~ 10 11 M© have essentially a 
spheroidal shape without major disk contribution whereas the lower 
mass systems show a large range in axial ratios. 

The average [a/Fe] gradients are consistent with zero for the 
slow rotators while a significant fraction of the low-mass, fast rota- 
tors show mildly positive gradients which we interpret as signs of 
the contributions from the central, more metal rich and abundance- 
ratio depressed regions with disk-like kinematics. The mean [a/Fe] 
gradients of the fast and slow rotators are 0.060 ± 0.012 and 
0.023 ± 0.016, respectively (the mean was determined with an out- 
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Figure 15. Metallicity gradients as derived from radial averages along 
isopho tes vers us f ex tra, the fraction of extra light as determined by 
iHopkins et alj (2009a) for galaxies in common. The gre y shaded area 
shows the approximate region occupied by the simulations of Hop kins et alj 
(2009a, their Figure 29). Fast and slow rotators are indicated by blue trian- 
gles and red circles, respectively. Cluster and field galaxies are shown as 
filled and open symbols. NGC4550, classified as slow rotator due to two 
counter-rotating disks, is shown as orange filled symbol. 



lier resistant bi-weight method and the errors reflect errors on the 
mean values). 

In the recen t simulations of major mergers of gas-rich disks 
bv lHopkins et alj (2009a). a large range of metallicity gradients is 
found. They emphasize the importance of the degree of dissipation 
in the central, merger-induced starburst for the resulting metallic- 
ity gradient. Their predictions (shown as grey shaded regions in 
Figures 

[13] and Ql broadly agree with our observations. This in- 
cludes the overall downturn in metallicity gradient with increasing 
mass. For galaxies in common we specifically explore the correla- 
tion between our observed metallicity gradients with their observed 
measurements of extra light f cx tra and compare this also with their 
predictions from simulations (grey shaded area) in Figure[l5] While 
the predictions broadly agree with the observations, a firm conclu- 
sion is not possible with this limited number of galaxies. 



4.3 Scaling relations of integrated stellar population 
estimates 

After the presentation of the stellar population maps (Section [4. It 
and the average gradients (Section l4~2t we now discuss more global 
SSP-equivalent stellar population parameters. Our integrated mea- 
surements were derived by averaging the (luminosity-weighted) 
spectra within a circular aperture with a radius of one R c and re- 
measuring the line strength indices. Small aperture corrections, as 
described in Paper VI, had to be applied for 28 galaxies. From the 
resulting line strength values (Table IaTI of the Appendix), we de- 
rived the stellar population estimates as described in Section \T2\ 
and presented in Table Q] We emphasize that these 'global' esti- 
mates of SSP-equivalent stellar population parameters derived from 



data covering essentially one effective radius (half the light) for 
nearby galaxies are much more representative of the galaxies as a 
whole compar ed to the typically central (often RJ&) values used in 
the p ast (e.g.. iTrager et al.ll2000l ; lKuntschnei]|2000l ; iThomas et al.1 
l2005t) . 

In Figure[T6]we show the SSP equivalent age, metallicity and 
[a/Fe] versus loga c , the dynamical mass Afdyn a nd log, A r re- 
lations for the one R c aperture. In agreement with IThomas et all 
(2005) we find a correlation of age with velocity dispersion a e . 
This correlation is driven by the very young, low-mass galaxies 
(<7 C < 100 km s -1 ) showing widespread signs of young stellar 
populations and the substantial fraction of fast rotators in the in- 
termediate mass range which are affected to various degrees by 
the contribution from younger stars linked to disk-like kinematics. 
Since the current sample is not complete, particularly at the low 
mass end, it is difficult to assess how many low-mass, old systems 
exist. Perhaps more importantly, the age versus dynamical mass 
plot shows that a strong contribution from young stellar popula- 
tions is confined to the low mass end (< 2 x 10 10 Mq). With in- 
creasing mass the relative fraction of 'young' galaxies decreases 
because any recent star-formation contributes relatively less to the 
total mass of the system. The slow rotators appear to be the oldest 
galaxies in the sample with a weak trend of decreasing \r and in- 
creasing SSP-equivalent age. Confirmation of this trend awaits the 
analysis of a larger more complete sample of slow rotators. How- 
ever, there is also a substantial number of fast rotators for which 
the last star-formation event was either long ago or very minor in 
mass fraction so that they show old SSP-equivalent ages. 



In agreem ent wi th previous 



Kuntschner et al. I2002I; If homas et alj 



2006; Collobert et al. 2006; Thom as et al 



inves tig ations (e.g., 



20051; iBernardi et ail 



2010) there is mild 



evidence for early-type galaxies residing in low density regions 
to have slightly younger ages on average compared to the cluster 
environment. For the SAURON sample, we find mean ages of 
10.8 ± 0.8 and 9.7 ± 0.9 Gyr, for cluster and field galaxies, 
respectively (the errors are given as errors on the mean). 

In summary, massive slow rotators have typically old SSP- 
equivalent ages while low-mass fast rotators make up the most 
prominent examples of young SSP-equivalent ages. In the inter- 
mediate mass range, galaxies show a range of SSP-equivalent ages 
of greater than 5 Gyr. The contributions from younger stellar pop- 
ulations are linked to secondary star-formation in disk- or ring-like 
structures (see Section |4.1.1t . 

Despite the various, distinct components and structures we 
find in the stellar population maps and discuss in Section |4~T1 the 
global, integrated metallicity estimates over one effective radius 
show a significant, positive correlation with log a c . A linear fit, ex- 
cluding NGC 3032, gives (see Figure [T6t: 

[Z/H](7? ) = 0.32(±0.07) loga c - 0.82(±0.15) . 

The metallicity - log a e relation is in agreement with pre vious in- 
vestigations (e.g., iKuntschnen'boOOl ; IThomas et ai]|2005l) . On av- 
erage we find smalle r metallicities compared to the study of e.g., 
Thom as et al. I d2005h , which is caused by the larger aperture we 
use in our analysis (one R c versus R e /10). The most extreme case 
of a post star-burst galaxy in our sample, NGC 3032, deviates from 
the relation. However, we ascribe this to the extreme luminosity 
weighting of young and old stars in this galaxy and therefore con- 
sider the SSP-equivalent estimate of the metallicity as not mean- 
ingful. There is no evidence for an environmental dependence of 
the metallicity-CT c relation in our sample. 

The correlation of our global (one R a ) metallicity estimates 
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Figure 16. One effective radius, integrated measurements of SSP-equivalent age, [Z/H] and [a/Fe] versus log a c dynamical mass M^ yn and log \r relations. 
Fast and slow rotators are indicated by blue triangles and red circles, respectively. Cluster and field galaxies are shown as filled and open symbols. NGC4550, 
classified as sl ow rotator due t o two counter-rotating disks, is shown as orange filled symbol. The dot-dashed lines in the Age and [a/Fe] versus log cr c relations 
are taken from lThomas et alj J2005l) . whereas solid lines are linear fits to our data. 



with dynamical mass Mdyn is significant at the 3a level but shows 
a significant scatter in metallicity at a given mass (rms = 0.08 dex): 

[Z/H] = 0.06(±0.02) log Mdyn - 0.75(±0.25) . 

iTrager e t al.1 feood) showed in a study of 50 local elliptical galax- 
ies that in the four-dimensional space of metallicity, log age, log a 
and [a/Fe] only a two-dimensional plane is occupied. Furthermore, 
velocity dispersion seems to be the only structural parameter that 
determines the stellar population of a galaxy. One of the projections 
from this plane is the remarkably tight correlation between [a/Fe] 
and log a. We show our [a/Fe] estimates as function of central ve- 
locity dispersion in the bottom left panel of Figure [T6] In good 
agreement with the literature, we find a rather tight relation in the 
sense that more massive galaxies exhibit l arger a/Fe ratios (see also 
j0rgensenlll999l;lKuntschner et alJlioOlt IProctor & Sansomll2002l ; 
Thomas et al.ll2005l l201fjT We do not detect a significant differ- 
ence between low and high density environments. The best fitting 
relation for the full sample is: 

[a/Fe] = 0.35(±0.06) loga c - 0.53(±0.13) . 

The best fitting relation of [a/Fe] with dynamical mass is: 

[a/Fe] = 0.09(±0.02) log M dyn - 0.70(±0.20) . 

What drives the [a/Fe] - log<r c (or mass) relation? Taking our 
results from Section 14.11 into account we find that the fast rota- 
tor galaxies with intermediate velocity dispersions (cr c = 100 — 
160 km s" 1 ) are the ones with contributions of stellar populations 
formed in a component with disk-like kinematics and more solar 
[a/Fe] values. The luminosity-weighted averages of the disk-like 



component with the main body of the galaxy results in interme- 
diate values of [a/Fe]. The group of galaxies with SSP-equivalent 
ages of < 4 Gyr have the strongest relative contribution of a young 
disk and SSP-equivalent estimates of [a/Fe] ~ 0.1 and thus deter- 
mine the low-mass point of the relation. We conclude that in our 
sample there is evidence that the young stars with more solar-like 
[a/Fe] ratios, created in fast-rotating disk-like components in low 
and intermediate mass galaxies, reduce the global [a/Fe] and thus 
significantly contribute to the apparent [a/Fe] - logo" e relation. 
However, we also find at least one low mass galaxy (NGC4458; 
a c < 100 km s _1 ), which has a [a/Fe] value comparable with the 
most massive systems in our sample. 

For the slow rotators, there is evidence for a trend of increas- 
ing [a/Fe] with decreasing \r. Overall, however, the metallicity 
and [a/Fe] estimates of the slow rotators appear surprisingly inho- 
mogeneous. Although the current sample of slow rotators is lim- 
ited, the significant spread in the detailed stellar population proper- 
ties may indicate that there could be more than one formation path 
resulting in slow rotators. 



5 CONCLUDING REMARKS 

We caution the reader at this point that the SAURON sample is a 
representative sample of limited size and not complete; therefore 
results for the population of early-type galaxies as a whole await 
confirmation with larger samples. 

Nevertheless, the SSP-equivalent stellar population maps pre- 
sented in this paper show that the 48 early-type galaxies in the 
SAURON sample display a significant and varied structure in their 
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stellar populations. Up to 50 per cent of the sample galaxies show 
signs of a contribution from a young stellar population. The young 
stars, dominating the light, are very apparent in the low mass 
(< 3 x 1O 1O M0) post-starburst systems, showing signs of resid- 
ual star-formation in the NUV and mid-IR and even optical gas 
emission lines consistent with star-formation in the most extreme 
case (NGC 3032). For these systems we observe SSP-equivalent 
ages of ^3 Gyr over the full field-of-view, typically covering one 
effective radius R a . Spatially well localized cases of young stars 
in circumnuclear disks and rings are found in eight intermedi- 
ate mass systems. The latter star-formation can be linked to thin, 
dusty disks/rings also seen by star-formation signatures in the NUV 
or mid-IR and optical gas emission lines consistent with star- 
formation for two galaxies (NGC 4459, NGC 4526). 

About 30 per cent of the sample galaxies show intermedi- 
ate SSP-equivalent ages (3-9 Gyr) where ongoing star-formation 
is rarely detected by NUV or mid-IR signatures, suggesting that 
the star-formation process was concluded more than ~1 Gyr ago. 
Common to most of these galaxies is the ubiquitous presence of 
ordered stellar motion, classifying them as fast-rotators. We show 
that the flattened components with disk-like kinematics identified 
in all fast rotators (Paper XII) are connected to regions of distinct 
stellar populations. This is clearly seen in the examples of the very 
young, still star-forming circumnuclear disks and rings with in- 
creased metallicity preferentially found in intermediate-mass fast 
rotators (e.g., NGC 4526). However, signs of increased metallicity 
and mildly depressed [a/Fe] ratios compared with the main body 
of the galaxy are seen in many fast rotators extending to apparently 
old galaxies with extended disk-like kinematics (e.g., NGC 4660). 
Additionally, intermediate mass fast-rotators show often evidence 
for nuclear stellar disks with increased metallicity. A galaxy for- 
mation picture emerges where essentially all fast rotators have ex- 
perienced, often multiple, periods of secondary star-formation in 
kinematically disk-like components on top of an older, spheroidal 
component. 

In contrast, the slow rotators, often harboring kinematically 
decoupled components in their central regions, generally show no 
stellar population signatures over and above the well known metal- 
licity gradients in early-type galaxies and are largely consistent 
with old (^10 Gyr) stellar populations. 

Using our radially averaged stellar population gradients we 
find that low mass fast rotators form a sequence of increasing metal- 
licity gr adient with in c reasin g mass, consistent with the observa- 
tions of ISpolaor et al] (§009). For more massive systems (above 
~3.5 xlO lo M ) there is an overall downturn such that metallic- 
ity gradients become shallower with increased scatter at a given 
mass leading to the most massive systems being slow-rotators with 
relatively shallow metallicity gradients. The shallower metallicity 
gradients and increased scatter could be a consequence of the com- 
petition between different star-formation and assembly scenarios 
following a general trend of diminishing gas fractions and more 
equal mass mergers with increasing mass, leading to the most mas- 
sive systems being devoid of ordered motion and any signs of recent 
star-formation. 

We also use our observations to compute more global stel- 
lar population parameters integrated over circular apertures of one- 
eighth and one effective radius, and compare the resulting relations 
of stellar population parameters versus velocity dispersion with 
previous results. In agreement with previous investigations we find 
significant metallicity - a and [a/Fe] - a relations. The metallicity 
- dynamical mass relation is shallow and shows a large scatter at a 
given mass possibly reflecting the many different formation and as- 



sembly histories which can lead to an early-type galaxy of a given 
mass. We argue that the [a/Fe] - a relation is largely driven by an 
increasing relative contribution from metal rich and mildly [a/Fe] 
depressed disk-like components in galaxies with decreasing mass. 
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APPENDIX A: REVISED LINE STRENGTH 
MEASUREMENTS FOR CIRCULAR APERTURES OF 

.Re/8 AND R E 

We present revised line strength measurements on the Lick/IDS 
system over circular apertures of i? c /8 and R c in Table lAll This ta- 
ble supersedes the data presented in Table 5 of Paper VI. The details 
of the revised data-reduction are given in Section[2] We further note 
that the effective radius values were revised using a R 1/n growth 
curve analysis of our wide-field MDM (1.3m) imaging survey of 
the early-type galaxy sample (Falcon-Barroso in preparation). 
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Table Al. List of revised line strength measurements within a circular aperture of Re/8 and Rc. 
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A 7A 
4. /U 


O 1 1 
Z. 1 1 


7 17 
Zl / 


1/1/1 
1 .44 


A SA 
4.31 


A 1 7 
4.1 / 


\jrjr A 17zl 
INLrL 43 /4 


7A 7 
/U.Z 


000 

zyz 


1 IS 
1.35 


S Aft 
3.U5 


A ^0 
4.3Z 


OA 
Z.Z4 


0A1 
ZOl 


1 AA 
1 .4U 


A A1 
4.D1 


A A7 
4.U / 


Mpr /I ICO 
INLrL 43 oZ 


O/l /I 

y4.4 


1 Q7 
15 / 


O 1 1 
Z. 1 1 


^ AA 
3.U0 


1 AS 
3.43 


O 1 1 

Z. 1 1 


1 7ft 
1 /5 


1 OQ 

i .yy 


A AA 
4.0U 


1 07 
3.Z / 


Mfp /1 107 
INLrL 435 / 


1 1 fi 
1 l.U 


1 A A 

1UU 


1 AO 

1 .oz 


S A7 

3.U / 


i o< 
3.y3 


O OO 

Z.ZZ 


QQ 

yo 


1 SA 

1.34 


A SA 
4.34 


1 AO 

3.oy 


Mf f" 1 /I /I ^ Q 
INLrL 4435 


ion 

iy.y 


1 AA 

1U0 


1 A1 

1.01 


A AA 
4.0U 


1 QO 

3.5Z 


1 QO 

i.yz 


Q1 
53 


1 A1 

1.01 


A r\S 
4.U3 


i ii 

3.33 


virir AA^Q 
INLrL 443 y 


A1 A 
4 l.U 


1 7G 
1/5 


1 A 
Z. 14 


^ 0^ 

3.Z3 


1 ft'N 

3.53 


OA 

z.zo 


133 


1 SA 
1 .54 


A AA 
4.DU 


1 17 
3.3 / 


Mpr A A 71 
IN LrL 44 / 3 


OA ft 
Z0.5 


1 OO 

iyz 


1 ACi 

i .4y 


S 17 
3.3 / 


A ^1 

4.31 


O 11 
Z.3 1 


1 QA 
150 


1 CA 
1 .3U 


A 7Q 
4. /5 


1 OQ 

3.y5 


Xinr AA77 
IN LrL 44 / / 


AA ^ 
40.3 


1 7A 
1 l\J 


1 AO 
1 .OZ 


s, AA 
3.U0 


A 1A 
4.14 


00 
Z.ZZ 


1 A7 
14 / 


1.33 


A A7 
4.0 / 


1 S1 
3.51 


Mfr A A OA 
INLrL 44 50 


1 AA O 

1U0.Z 


111 

311 


1 1 < 
1.13 


3.13 


C 1 1 

3.11 


0< 

Z.Z3 


OAQ 
Z05 


11^ 
1 .13 


A Sft 
4.3y 


A Sd 

4.3y 


INLrL 43 ZO 


1^ 7 
33. / 


010 
Z3Z 


1 ft/1 
1 .54 


^ ^ 1 

3.3 1 


A 11 
4.3 1 


10 
Z.3Z 


01 A 
Z14 


1 ^S 
1.35 


A Q1 
4.y 1 


A 1 A 
4. 10 


Mf /- 1 /I ^ A A 
INLrL 4340 


00 A 
ZZ.U 


00A 
ZZU 


1 SA 
1.34 


S 1/1 
3.34 


A AA 
4.0U 


OQ 
Z.Z5 


1 QO 

ioy 


1 SA 
1 .34 


A <A 
4.30 


1 Q7 
3.5 / 


INLrL 43 3 U 


1 1 A 
11.0 


55 


1 1 
Z. 1 1 


A ^A 
4.34 


1 1 7 
3.1/ 


AA 


1 A1 
1U3 


1 QO 

i.yz 


A 1^ 
4.33 


QA 

z.yo 


Mrr /l ^^o 
INLrL. 43 3Z 


11 O 

33. y 


OAQ 
Z05 


1 SS 
1.33 


S A^ 
3.03 


A ftl 
4.53 


OO 

z.zy 


Oil 

Z33 


1 AS 
1.43 


S 1 A 

3.14 


A /1A 
4.4U 


Mfp /I ^ A/1 
INLrL 43 04 


1 Q 1 

iy.3 


1 71 
1/3 


1 ^o 
1.3Z 


^ Q7 

3.0 / 


/I 77 
4. / / 


O /lO 
Z.4Z 


1 <A 

13U 


1 ^7 
1 .3 / 


A Q/1 
4.54 


/I AA 
4.UU 


NGC4570 


12.8 


197 


1.45 


5.80 


4.67 


2.33 


167 


1.45 


4.75 


3.97 


NGC4621 


46.0 


225 


1.40 


5.29 


4.68 


2.27 


200 


1.44 


4.67 


4.11 


NGC4660 


ll.5 


221 


1.43 


5.56 


4.76 


2.30 


181 


1.47 


4.71 


4.02 


NGC5198 


18.0 


206 


1.46 


5.44 


4.73 




173 


1.55 


4.72 


3.97 


NGC 5308 


9.9 


244 


1.46 


5.28 


4.46 




201 


1.46 


4.83 


4.16 


NGC5813 


55.9 


226 


1.49 


5.00 


4.47 


2.12 


210 


1.49 


4.63 


4.03 


NGC 5831 


29.2 


168 


1.72 


5.55 


4.13 


2.38 


148 


1.81 


4.83 


3.34 


NGC 5838 


20.6 


283 


1.59 


5.68 


4.60 


2.38 


232 


1.57 


4.96 


4.17 


NGC 5845 


4.3 


269 


1.51 


5.62 


4.64 


2.30 


237 


1.51 


5.23 


4.36 


NGC 5846 


76.8 


236 


1.33 


5.37 


4.80 




213 


1.30 


4.83 


4.34 


NGC 5982 


24.9 


260 


1.62 


5.74 


4.41 




223 


1.59 


5.11 


3.95 


NGC 7332 


9.2 


137 


2.27 


5.87 


3.67 


2.28 


125 


2.14 


4.85 


3.19 


NGC 7457 


33.2 


70 


2.37 


4.98 


2.92 


2.09 


75 


2.25 


4.56 


2.80 



Notes: This table supersedes the values published in Table 5 of Paper VI. (1) NGC number. (2) Effective (half-light) radius R c measured with a ij 1 /" growth 
curve analysis from our wide-field MDM (1.3m) imaging survey of the early-type galaxy sample (Falcon-Barroso in preparation). (3) Velocity dispersion 
of the luminosity weighted spectrum within a circular aperture of R c /S. Only the first two moments, v and a are used in the fit to the spectrum. (4) - (7) 
Fully corrected line strength index measurements on the Lick/IDS system of the luminosity weighted spectrum within R c /S for the H/3, Fe5015, Mg b and 
Fe5270g indices. Due to limited field coverage we cannot determine the Fe5270s indices for 1 1 galaxies. (8) Velocity dispersion of the luminosity weighted 
spectrum within a circular aperture of one R c . Only the first two moments, v and a are used in the fit to the spectrum. (9) - (1 1) Fully corrected line strength 
index measurements on the Lick/IDS system of the luminosity weighted spectrum within R c for the H/3, Fe5015, and Mg b indices. For galaxies with less 
than one R c coverage we applied the aperture corrections given in Paper VI. Formal errors of the line strength indices are below 0. 1 A. We note, however, 
that there are systematic errors which we estimate to be of the order of 0.06, 0. 15, 0.08, 0.06 A for the H/3, Fe5015, Mg b, and Fe5270 s indices, respectively. 
For the velocity dispersion we adopt an error of 5 per cent. 



